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CONSERVATION LAWS OF LINEAR, HOMOGENEOUS SYSTEMS 

By J. R. Williams 
Electronics Research Center 


SUMMARY 


Conditions for a linear, uniform, homogeneous system to have 
nonsingular quadratic invariants, formed in the usual manner from 
constant metrics, are expressed in terms of properties of a 
Kronecker sum in the operator. Useful tests for the existence of 
such invariants are described. It is shown that system definition 
in terms of these invariants is not only restricted in its appli- 
cation but inconclusive. The interpretation of quadratic invari- 
ants, formed from metrics that depend on the space variable, as 
conservation laws is supported both formally and by physical argu- 
ments for the nonuniform system. A property of the Kronecker pro- 
duct and sum, that the Kronecker product of the matricants of two 
system operators is the matricant of the Kronecker sum of the opera- 
tors, is then utilized to describe a useful procedure for obtain- 
ing the metrics, including the constant metrics. This property 
has application also in group-representation theory, and is used 
to obtain a representation of the abstract group onto a group of 
matrices that display the metrics defining the group. 

The operator for a linear, lumped system is treated within 
the framework of the theory. Techniques developed for uniform 
distributed systems are found to be appropriate, and a group- 
theoretical interpretation is given. 


I . INTRODUCTION 

In this report we present a theory of the conservation laws , 
expressed as quadratic invariants, obeyed by a linear, homogeneous 
system. The theory applies to a distributed system that can be 
described by a matrix operator, over the field of complex numbers, 
relating a state vector and its derivative in a single space vari- 
able; and to a lumped- system operator relating a state vector at 
two generalized ports. The distributed-system operator is allowed 
to be a function of the space variable under a weak assumption of 
bounded differentiability, thereby admitting nonuniformity in the 
system due to any physical mechanism. The nonuniformity could be 
due, for example, to nonuniform boundaries in a wave-guiding sys- 
tem, material inhomogeneity, or nonuniform coupling in a single- 
frequency or parametric device. New insight into the quadratic 
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invariants of uniform systems is gained when these are considered 
in the context of the more general theory, and extended techni- 
ques for their study are developed. Operators for distributed 
and lumped systems are included in a unified theoretical treat- 
ment. We assume that the system has a finite number of degrees of 
freedom, but the fundamental results could be generalized. 

Important results concerning the nonsingular quadratic invari- 
ants of linear, uniform, homogeneous systems, and of linear lumped 
systems, have been given by Pease (refs. 1, 2, and 3). A space 
is considered to be defined by a metric, which specifies how the 
length of a state vector is to be measured, and a conservation 
law to be determined by each metric for which the operator becomes 
rotational. Under certain conditions on the eigenvalues and 
structure of the operator, the system will obey a set of conserva- 
tion laws determined by linearly independent nonsingular metrics, 
equal in number at least to the number of degrees of freedom of 
the system. The proof of this result provides a formal procedure 
for obtaining the metrics, a procedure that requires a determina- 
tion of the eigenvalues, eigenvectors, and generalized eigen- 
vectors of the operator. The metrics for which the operator be- 
comes rotational can also be obtained by writing out and solving 
the governing equations. A set of linear, homogeneous equations 
is obtained, where the coefficient matrix is of squared order 
compared with the operator and the metrics. Either approach can 
be tedious if the operator matrix is not of small order. 

Conversely, given two metrics of suitable form, a range of 
distributed-system operators that become rotational for either 
of these metrics can be determined. This range will in general 
be narrowed if more metrics are specified. The distributed sys- 
tem can be represented by a matricant — obtained by integrating 
the differential equations describing the system — which relates 
the state vector as a function of position to a boundary value. 

The specified metrics define a group for the abstract operator, 
the matricants determined by the metrics are a representation of 
the abstract group onto a Lie group, and the system operators are 
members of a corresponding Lie algebra. From this point of view, 
systems are defined by their conservation laws. We shall discover, 
however, that in general an operator is not determined uniquely 
by specifying all of the metrics for which it is rotational; in 
fact, the simplest class of operators for conservative systems 
will be determined only within a space of dimensionality equal 
to the number of degrees of freedom of the system. Most opera- 
tors do not become rotational for any metric. 

The coefficient matrix, obtained by writing out the equa- 
tions governing the metrics, while of squared order, signifi- 
cantly takes the form of a Kronecker sum. This has some practical 
as well as important theoretical consequences. First of all, it 
will permit us to develop a formal parallel to the theory developed 
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for the operator, in terms of properties of the coefficient matrix. 
Existence criteria for nonsingular solutions, as well as other in- 
formation about the operator and the metrics , are expressed in 
terms of more accessible properties of the coefficient matrix; 
and equivalent procedures for obtaining this information for an 
operator in the restricted class considered can be compared with 
those available for testing the operator directly. 

With the theory reformulated in this way, we shall be able 
to complement the previous results by introducing metrics that 
vary with the space variable, and to confirm that these are re- 
quired for a definitive theory of this and the larger class of 
operators . Quadratic invariants formed from such metrics , in 
general, can meaningfully be interpreted as conservation laws, 
and uniform systems can be considered a kind of "singularity" 
for which this interpretation breaks down. These and the constant 
metrics can be obtained in a way that indicates that they, as 
well as the matricants , have overt group-theoretical properties. 

Not only is a useful procedure for obtaining the metrics revealed, 
but a representation of the abstract group explicitly in terms 
of the metrics is indicated. The Kronecker sum, formed from the 
distributed-system operator, will then be recognized as the opera- 
tor in the corresponding Lie algebra. Such a representation will 
be seen to be what is required for definitive system definition 
for the uniform system as well, which can be viewed as a limiting 
case when interpreting the invariants as conservation laws . 

We shall also discuss briefly another class of metrics that 
depend on the space variable, and that determine ranges of non- 
uniform systems in a manner very similar to that in which the 
constant metrics determine ranges of uniform and nonuniform sys- 
tems. While there is some basis for interpreting quadratic forms 
obtained from these as conservation laws, and while similar tech- 
niques for synthesis from these metrics can be employed, they have 
the same limitations for system definition as do the constant 
metrics . 

Operators for lumped systems, giving as they do "black box" 
descriptions, are less definitive in circumscribing the kinds' of 
systems they might describe. Their invariants have been obtained 
in a manner parallel to that for distributed systems, but this as 
well as the synthesis problem can in fact be reduced to the prob- 
lem for uniform distributed systems. From the point of view of 
group-representation theory, this reduction in the range of sys- 
tems that need be considered can be viewed as resulting from a 
homomorphic mapping of the group of abstract operators for the 
more general nonuniform system onto the full linear group. 

The material presented in this paper, like the results ob- 
tained previously, has broad applicability. So that it may be 
most generally useful, the presentation is abstract. No attempt 
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is made to discuss sophisticated physical applications; but guide- 
lines for the physical interpretation of results are indicated, 
and sufficient illustrative material is included to make the theo- 
retical points clear. 


II. BASIC CONCEPTS 

In this section we summarize, for review and reference, some 
results concerning the quadratic invariants of linear, uniform, 
homogeneous systems . The analytical techniques by which these 
results have been obtained are outlined, for comparison with the 
approach in Section III, and several topics are pursued in more 
depth. Finally, we call attention to some tests for such invari- 
ants that are equivalent to applying directly the existence 
criteria that emerge in the formal development of the theory . 

The converse problem to the central problem discussed in this 
section is perhaps of even more fundamental importance . We re- 
view current concepts in dealing with this problem in Sections IV 
and VII, at more logical points in the development of this paper. 


The System 

We assume first a linear, distributed system having a finite 
number of degrees of freedom, such as can be described by a first- 
order differential equation of the form 

^ = - jRx (1) 


Here x is a column matrix, an n-dimensional state vector describ- 
ing the state of the system as a function of position, and JR is 
an n*n matrix, the system operator. The components of x^ are con- 
sidered to be functions of the single space variable z, and the 
time dependence has been removed. Such a system is called homo- 
geneous, for Eq. (1) is homogeneous in x^. 

As a simple example, on a voltage-current basis. 



the system operator for a uniform transmission line of character- 
istic impedance Zq , having an associated propagation constant 3/ 
is 
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0 


ez o\ 

] (3) 

V 

The form of Eq. (1) is not as restrictive as may at first 
appear. It is not restricted altogether to the description of 
systems operating at a single frequency. Furthermore, order can 
be traded for dimensionality; and accordingly all systems of lin- 
ear differential equations, of finite order and dimensionality, 
can be reduced to this form. 



Metrics and Quadratic Forms 


Before discussing the conservation laws governing such a sys- 
tem, we first introduce the concept of a quadratic form , or gen- 
eralized inner product. This is a scalar quantity defined by 


s (K) 


+ 

x Kx 


(4) 


where the superscript t denotes the complex conjugate of the trans- 
posed vector, and where K is a constant nxn matrix. 

In order to suggest the appropriateness of such forms for 
describing physical properties of a system, we note that if the 
components of x are voltage and current, as in Eq . (2), then with 



the power may be written 



We postpone further consideration of the physical signifi- 
cance of these forms — except to point out that K is taken to be 
constant so that the interpretation of s (K) will^not depend on 
position. 
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.j. Equation (4) is a generalization of the unitary inner product 
xjx, for which & is taken as the identity matrix 1^, and provides 
a more general definition of the square of the length of the vec- 
tor x. The matrix IK is called the metric of the space in which x 
is said to be embedded, and as such can be considered to define 
the space . 


Appropriate Metrics and Invariant Forms 

Of particular concern to us here will be those metrics K for 
which s (KJ is invariant under the operation jR. In other worcfs , 
for a given R, we are interested in those spaces, if any, in which 
x can be embedded, such that its length will not vary with z. 

Xs sociated with each such then is a conservation law, which pre- 
dicates a quadratic form that remains constant as x^ varies in 
accordance with Eq. (1) . We shall refer to these K as the metrics 
appropriate to JR. 

For s (K) to be invariant under the operation R, it is required, 
by Eqs . (l^and (4), that w 


ds _ d%^ 
dz dz 


Kx + x + K 


dx 

iS£S 

dz 


jx + (R + K-KR) 


x = 0 


(7) 


+ 

where R is the transpose conjugate of R^. Since Eq. (7) must hold 
for any x. , for to be appropriate to _R~we must have 

R + K^ = KR_ (8) 


If K is a nonsingular solution to Eq. (8) , R is said to be K- 
Hermitian . ~~ ~ 

Thus far a restriction has not in fact been imposed on the 
elements of R. that they may not be functions of z. On the assump- 
tion that the system is uniform, however, and R_ is thus constant, 
results have been derived concerning the existence of appropriate, 
nonsingular jK. Singular jK have not been considered important to 
the theory. 


Necessary Conditions for Appropriate, Nonsingular 

+ 

For Eq. (8) to hold for some nonsingular K, R and R must be 
related by a similarity transformation, , ~~ 

R + = KRK -1 (9) 

MM Ml M MOM 
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This imposes conditions on the eigenvalues and structure of R^. 

The n eigenvalues of may be distinct, or some may be re- 
peated. Associated with each distinct eigenvalue Bi is a linearly 
independent eigenvector x(i/l), satisfying 


Rx 


(i / 1) 




(i / 1) 


( 10 ) 


Associated with an eigenvalue repeated m^ times are nu linearly 
independent eigenvectors and generalized eigenvectors, forming 
one or more chains. We shall use different subscripts to index 
repeated eigenvalues when these are associated with independent 
eigenvectors or separate chains, but we shall denote the number 
of eigenvalues equal to any Bi by rtij_ (i.e., if Bi = 6 j , mj_ = m j ) . 
The chain of length , l<S.i<mj_, associated with the eigenvalue 
Bi, consists of an eigenvector , satisfying Eq . (10), and if &j_>l, 

of generalized eigenvectors x(i'k), k = 2 , 3 , . . . , , satisfying 


Rx (i ' k) = B.x (i ' k) + x (i ' k - 1) 

J. * VWV 


(ID 


( -j V ) , 

where x. v ' ' is said to be of rank k. A matrix is said to be 
semisimple if it has a complete set of n linearly independent 
eigenvectors — i.e., if the chains are all of length 1. It is 
necessarily semisimple if its eigenvalues are distinct. 

Now it is readily shown by Eqs . (10) and (11) that if two 

matrices A^ and are related by a similarity transformation, A_ = 

S“1bS , and if A has a chain x(ifk) k = 1,2,. with eigen- 

value Bi / then B has a chain Sx'' 1 ' K ' , k = l,2,...x.j, with eigen- 
value Pi* 

Thus necessary conditions for the existence of a nonsingular 
solution K to Eq. (9) are that R and have the same eigenvalues 
and structure. But the eigenvalues of ^R T are the complex con- 
jugates of the eigenvalues of ^R. Also, if R^ has a chain , 

k = 1 , 2 , . . . , &i , with eigenvalue Bi , R^ will have a chain / 

k = 1 , 2 , . . . , , with eigenvalue B* ^ where the (jw(i/k))+ are 

the reciprocal eigenvectors and generalized eigenvectors of R^ 
with eigenvalue Bi • The conditions then become that the eigen- 
values Bi of R either be real, 


Bi = B* 


(12) 


or occur in conjugate pairs. 


B • = B • 


(13) 
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such that chains of equal length are associated with 6^ and g ^ . 


Sufficiency of Conditions, Degeneracy, and Formal 
Solution — R^ Semisimple 


That these conditions on the eigenvalues and structure of R 
are also sufficient, and that there will then be at least n lin~ 
early independent, nonsingular appropriate to R^, can be shown 
by exhibiting the appropriate K T In so doing, a"~formal procedure 
for their determination is described. 

On the change of basis 

x = Sx' (14) 


where S is a constant, nonsingular n*n matrix, Eq. (1) becomes 


dz 


. _ i i 

-jR x 


(15) 


where 


' -1 
R = S RS 

»vv\ /vv» wn 


(16) 


The quadratic form, Eq. (4), becomes 


(K ) = x' 


'K x 

MV 


(17) 


where the metric on the new basis is 

K' = S + KS (18) 

/va rv*% iwk 

Since S* is constant, Jk' is constant. Accordingly, s' ( k' ) will be 
invariant under the operation r' if w ' 

AAA 

I 4. I 1 1 

R ~K = K R (19) 

AAA AM AAA W) 

If so, s (K) will be invariant under the operation R, since s (K) is 
independent of the basis. ^ /W1 

In particular, S can be chosen to be a modal matrix of R — 
a nonsingular matrix?* the columns of which consist of the n Tin- 
early independent eigenvectors and generalized eigenvectors of R 


8 



(the complete set of eigenvectors, if II is semisimple) , with each 
chain arranged in the order of increasing rank. Since jR is con- 
stant, S can be chosen constant. 

AAA 

Assume first that R is semisimple. Then under the transforma- 
tion of Eq . (16), if S^Ts a modal matrix, R will be diagonalized 
on the new basis. Thus 



where (3^/ appearing in the ii position, is the eigenvalue t that 
corresponds to the eigenvector x^(i / 1) appearing as the i th column 
of S. If the are real or occur in conjugate pairs, the condi- 
tion* on the structure is satisfied automatically for R semisimple, 
and r'+ will be obtained from R* by interchanging the positions 
of conjugate pairs of eigenvalues. This interchange of positions 
can be effected by a similarity transformation, 

R ,+ = P -1 r'p (21) 

AAA AAA /V/V\ AAA 

where P = P -1 is a permutation matrix having a 1 in the ii posi- 
tion if* 6^= or in the ij and ji positions if is complex 

and Bi = 6 j • Then Eq. (19) becomes 


r' (PK 1 ) = (Pk' ) r' (22) 

AAA />M AM AM am am 


Any diagonal matrix will commute with Jjl . Accordingly, there 
will be at least n linearly independent metrics K appropriate to 
R ; namely, 

Aaa 



(23) 


where a.,a 0 ,...,a are arbitrary constants. Furthermore, if R has 
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repeated eigenvalues , then for each 3j_ = 3 j two additional arbi- 
trary constants are introduced, in the ij and ji positions of PK' . 
The p constants, p>n, can clearly be chosen to obtain p linearly" 
independent, nonsingular . 

Finally, a set of p linearly independent, nonsingular 
appropriate to .R is obtained by transforming back to the original 
basis. Inverting Eq. (18), we obtain for each K 1 , 

flAA 


K = (S + ) 1 k’s 1 

/W\ /vA AAa 


(24) 


Since S and K are nonsingular, K is also; 

AAA AAA r AAA | 

dence of the K follows from that of the K . 


AAA 


AAA 


and the linear 


indepen- 


,R Not Semisimple 

More generally, if S^ is a modal matrix of R, R is transformed 
by Eq. (16) to a quasidiagonal matrix in the Jordan' canonical form. 


R 


AAA 


quasidiag (J^,^, 



(25) 


a partitioned matrix in which, corresponding to each chain of 
length associated with the eigenvalue 3j_/ there appears an £j_x 

Jordan block, 


J.= 


/ 


e i 




0 

I 

I 

I 

I 

0- 



(26) 


positioned in the same columns of r/ as the chain is positioned 
in S. By assumption, either 3-i is real, in which case 

AAV 


j . 

AAA 1 


J* 

AAA 1 


(27) 


or there exists an £ . x£ 

3- = 3*/ such that 1 
j l 


I 

Jordan block of R 


AAA 


with complex eigenvalue 
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(28) 


J. = J. 


,-i 


We define a f permutation matrix V* = P partitioned in the 
same manner as . Indexing the positions of the component sub- 
matrices by rows and columns of blocks, we put an identity 

block I 0 in the ii position of P (the position of J. in r' ) if 

A/v '<~ ^ A/W M*± AaA 

* 

Oh , or in the ij and ji positions if is complex and = 


Ji = 

J. , . 

AA*- 1 - 

/WAX 

J* . 

/WIJ 

Then 


— 1 I 1 * 

P R P = R 

AM /w l AW AM 


= quasidiag (J, ,J„,...) 


since the positions in R of each Jj and satisfying Eg, 

• ... 1 , , , , . A ' w% , , . 'W\l maJ J ^ a 

are interchanged by this permutation. 


Let Qi^ be a modal matrix of . 


Then 


(29) 

(28) 


t * 

j!q. = Q. J . 

/Wl .AA*! 

• * *f* 

since J. is the Jordan canonical form of J.. Defining 

A*»l -**" 1 ^ 


AWl 


(30) 


Q = quasidiag (Q, ,Q 0 , . . . ) 

/vv\ sw-L r** Z 

it follows that 


(31) 


or 


Thus 


i 4* i * — I i 

R Q = QR = QP R P 

AW AAA AA WA AAA JMA. AAA 


r' (QP 1 ) = (QP 1 ) r' 

w> MAU aw 


K = QP 

WA AM AA 


-1 


(32) 


(33) 


is a solution to Eq. (19). 

.j. Each of the eigenvectors and generalized eigenvectors of 
J. , which make up the columns of Q., is specified only within an 

Mil 
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arbitrary constant. The expressions for Q and , then, each 
contain n arbitrary constants. Equation T33) reduces to Eq.(23) 
for JR semisimple, since postmultip lying instead of premultiplying 
by simply permutes the arbitrary constants. 

f 

Again, if R has more than one Jordan block with t the same 
eigenvalue, there will be further degeneracy in the K, . The 
blocks in Q off the quasidiagonal are not in general square, how- 
ever, and lEhe manner in which additional arbitrary constants are 
introduced becomes more complicated. 

As before, the p constants, p>n, can be chosen to obtain p 
linearly independent, nonsingular ^ ; and these can be trans- 
formed back to the original basis by Eq. (24) . 


Lumped Systems 

Similar results are obtained for a lumped linear system, 
described by the equation 

x~ = Mx. (34) 

i\A* AA* AML 


where M is an n*n nonsingular matrix relating a state vector at 
two generalized ports of the system. We seek metrics such that 

x^Kx = xTkx, (35) 

which requires that 

M + KM = K (36) 

/WV /w% <VV\ /Wl 

If K is a nonsingular solution to Eq, (36) , M is said to be K- 

/W\ J AW 

Unitary. ‘ 

Necessary conditions for M to be .K-Unitary are that the 
eigenvalues A^ of either be of unit magnitude, 

A. A* = 1 (37) 

l l 


or occur in pairs satisfying the cross-conjugacy relation 


A .A 


l 



(38) 
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such that chains of equal length are associated with A ^ and A j . 
Again these conditions are sufficient to ensure the existenceof 
at least n linearly independent K for which M is K-Unitary. 

AA/\ /tvs. av\ J 


The solutions to Eq. (36) are obtained as before by first 
transforming to M , the Jordan canonical form of M. Then by a 
similar procedure the K appropriate to M are obtained in the 
form of Eq . (33). The positioning of the blocks in the permuta- 

tion matrix depends in the same manner on the self-con jugacy or 
cross-conjugacy of the eigenvalues in the associated Jordan blocks. 
The Qi are, in this case, modal matrices of the 


These results for a lumped system can be applied directly to 
the analysis of a uniform distributed system also, if Eq. (1) is 
first integrated. If is nxn and constant, Eq. (1) will have n 
linearly independent solutions for x^. (The elements of R^, in 
fact, need only be differentiable with bounded derivatives in the 
region of interest — this is known as the Lipschitz condition.) 

If any set of linearly independent solutions is arranged as the 
columns of an nxn matrix M(z), then M(z) satisfies 


dM 

dz 


- jRM 

Am /wa 


(39) 


The solutions to Eq. (39) are known as integral matrices of R^. 
Since the columns of M(z) are linearly independent, we can form 
its inverse at any value z^ of z. Thus we can define 


M(z,z.) = M ( z ) M hz-,) 

1 AAA AAA 1 


(40) 


As the postmultiplication of M(z) is by a constant nonsingular 
matrix, M^ZjZi) is still a nonsingular solution of Eq. (39). It 
is a particular integral matrix having the property 


M(z. , Z ) 

/VVN i 1 


I 

Aaa 


(41) 


M(z,zi) is known as the matricant of the distributed operator R. 

Tn view of Eqs. (39) and (41) , w 

x (z) = M (z, z. )x (z, ) (42) 

A*A AAA 1 AAA 1 

is a solution to Eq. (1) , expressing x^ in terms of a boundary 
value. Equation (42) provides a lumped representation of the sec- 
tion of the distributed system between z-^ and z. 
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By comparing the eigenvalues and structure of and I4(z,zp) , 
for R constant, we could demonstrate the equivalence of the condi- 
tions for R to be K-Hermetian and for M(z,z-i) to be K-Unitary, 
with respect to the same number p>n of linearly independent K 
(except that the degeneracy of JM may be greater at discrete values 
of z) . If R is semisimple, M(z,zi) will be also, with the same 
eigenvectors; the eigenvalues of R are related to the eigen- 
values A . of M (z , zi ) by AA ‘ 

1 AV\ 

-jg . (z-z, ) 

A i = e 1 (43) 

The equivalence of Eqs . (12) and (13) , the conditions on the Bj_, 

and Eqs. (37) and (38) , the conditions on the Aj_ , is clear, and 
only at discrete values of z can Aj_=A j if Bi^Bj . If R is not semi- 
simple, the generalized eigenvectors of and M ( z , z± ) will not be 
the same — in fact, the generalized eigenvectors of M(z,z]_) will 
depend on z; but the chain of R with eigenvalue Bi will be re- 

-j6 i (z-z 1 ) 

fleeted into a chain of ^l(z,z]_) with eigenvalue e 

We do not need to compare details of the theory to see that 
we will obtain in fact the same appropriate metrics K by either 
approach. For let K be a metric appropriate to M(z,Zt), satisfy- 
ing Eq. (36) for all z. Differentiating, we obtain 

(M f KM) = kM + M + K = 0 (44) 

Q Z AAA /V«A r*\A Cl Z AMAAA Aaa AWV QZ AAA. 

Substitution of Eq. (39) and its transpose conjugate gives 

jM + (R + K-KR) M = 0 (45) 

/VVN AA'V AAA AAA AM AAA. AV* 

Since is nonsingular, we can premultiply by - j (M^ ) 1 and post- 
multipTy by M“l, yielding Eq . (8), which is the condition for to 

be appropriate to jR. The steps are also reversible, except that 
Eq. (44) requires only that JYMKM be constant. But since we can 
then evaluate it at any value o?" z, by Eq. (41) this constant 
matrix must be K, which proves the converse. 

Aaa 

We note for future reference that the preceding result 
applies as well if R^ is not constant, for any constant solutions 
to Eq. (8) , or to Eq". (36) if is the matricant. 

Thus any linear, distributed system, described by a system 
operator R satisfying the Lipschitz condition, has an equivalent 
lumped representation described by the matricant M(z,z^) of ^R, 
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the particular solution of Eq. (39) satisfying the boundary condi- 
tion, Eq. (41) . Only in special cases (including II constant) are 
there exact procedures for obtaining this solution. But if the 
matricant representation is given, it is simple to obtain R. In- 
verting Eq. (39) , 


R 


= 3 4m’ 1 

J dz A~ 


(46) 


The operators and jM(z,Zn) have the same appropriate metrics, 
constant satisfying* both Eqs . (8) and (36) . 

We also point out that a lumped system — for constant, non- 
singular M, not a matricant — can be represented equivalently by 
a uniforrrHdistributed system. We write 

M = SM's -1 (47) 

rS*\ AAA 

where S is a modal matrix of M, and 

mJ = quasidiag ( • • • • ) (48) 

is its Jordan canonical form. The £-j_x&^ Jordan blocks are given 
by Eq. (26), with Xj_ replacing 6 . For each Xj_ we define a Bi 
(which may be complex) as the solution to 


-jB.L 

e = X i , — 7T < ReB.^ it (49) 


where L is an arbitrarily chosen equivalent length, 
blocks 


/N 


id 


(Z) 



We define 


( — j z ) 



u ± -l) ! 
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These have the property 


j 


A 


dJ. (z) 

/vnl 

dz 



u . 


(z) 


J.(*> 


(51) 


where is the Jordan block of Eq. (26) ^ with eigenvalue 3i . 

Also, Jjl is the Jordan canonical form of J^ (L) , and J^(0) = I.^ • 
We then define i 


J • (z) 

/W»l 


A *i A A 

TT^J . (z ) T . 


(52) 


A A 

where T-: is a modal matrix of J-j (L) , and J-: (z) has the properties 

AAA- 1 - AM 1 - 


J ■ (L) 

amI 


J. 

aaaI 


J, (0) 

AAA _L 



(53) 


We form 


T = quasidiag (T, ,T„ , . . . ) 

/W\ aaa_L yv*Z 

If A A 

M (z) = quasidiag (J, (z) , J 0 (z) , . . . ) (54) 

AAA. AA*_L A/'/*/! 

m' (z) = quasidiag (J, (z) , J„ (z) , . • . ) 

rwx AaaZ 


It then follows that 

M ( z , 0 ) = SM'(z)S -1 = ST _1 m" ( z)TS -1 (55) 

/y/VN /Ws A*\ * /ywv AAA 

is the matricant of a uniform distributed system, with operator 


R = 

rtAA 



(ST 1 

/v** 


)r' (ST 1 ) 1 

AAA. AAA />M 


(56) 


where 


r' = quasidiag (jJ^ , ji^^ , . . • ) 

Mn .wil <wi 2. 

and satisfies the boundary conditions 


( 57 ) 
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( 58 ) 


M (L , 0 ) 

AW 

M ( 0 , 0 ) 

AM 


M 

I 

am* 


In general, S is a modal matrix of ^(Z/O) and ^ m' (z) is its 
Jordan canonical form only at z = L, but ST“1 is a modal matrix 

_36 i z 

of R. If M is semisimple, J^(z) = e . Then M^ (z) is the 

JorcPan canonical form of M(z,0) and S is a modal matrix of both 
M.(z,0) and 1^. In this case R is determined quite simply as R^ = 
< SR , S~1, witfT the 8i obtainecT^f rom Eq. (49). 

fW> /Wv /V* - 1 - 

The significance of this result is that we shall be able, 
without loss of generality, to develop the subsequent theory with 
a focus on distributed systems. Whatever the physical basis for 
the transformations they describe, lumped systems can be treated 
by techniques developed for uniform distributed systems. At this 
point we can observe that the metrics appropriate to M will be 
the constant solutions to Eq . (36) for the matricant ^f Eq. (55), 

and may therefore be obtained as the solutions to Eq. (8) for R^ 
in Eq. (56). That these will be all the appropriate is ensured 
by restricting Re8i in each case to the interval indicated in 
Eq . (49) — in fact, we need only require that 8i = 8* when A^Aj=l 


Hermitian Metrics 


Interest has centered on real invariant quadratic forms as 
representing real physical quantities conserved by the system. 

For s (Kj to be real, it is required that be Hermitian — i.e., 

K = K+ . But the K determined by the procedure described will 

AAO /WV AArt ... I , __ 

not in general be Hermitian. (Hermitian K do not imply Hermit- 

• T, \ 

ian K. ) 

nrr\ 

Hermitian K for which R is K-Hermitian or M is K-Unitary 
can, however, be formed from any set of appropriate, nonsingular 
K. If is appropriate to the operator, is also appropriate, 
aTs can be seen by taking the transpose conjugate of both sides of 
Eq. (8) or Eq. (36). Moreover, the appropriate define a space 
since Eqs . (8) and (36) are linear, linear combinations of appro- 
priate metrics are appropriate also. Thus from the given set 
we can form appropriate metrics 


3 a -; 


“Da 


K 

1 


= e 


K. + 

Mil 


Wl 



(59) 
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which will be Hermitian for and b-^ real. Furthermore, real a^ 
and b-i can always be chosen so that K and K, are nonsingular. 

1 X 

It may be shown that, of the 2p Hermitian metrics obtained 
by Eq. (59) from a maximal set of p linearly independent, appro- 
priate K^, exactly p will be linearly independent. Given a set 
of p linearly independent, Hermitian metrics appropriate to 

an operator describing a physical system, however, it does not 
follow that the physical significance of each of the associated 
quadratic forms will be apparent. But the span a subspace 

of the space spanned by the . Linear combinations ^2 c^K^^ will 

be Hermitian for c^ real, will include all appropriate, Hermitian 
metrics, and may generate quadratic forms having clear physical 
interest. The physical interpretation of the invariant forms of 
a system presents a separate problem in every individual case, 
beyond the determination of a linearly independent set of mathe- 
matically appropriate forms. 


By way of example, the system operator for a uniform trans- 
mission line, given by Eq. (3), has two distinct eigenvalues +0, 
which, assuming a lossless line, will be real. We therefore 
expect that there will be two linearly independent, nonsingular 
metrics K appropriate to R? and we find as an Hermetian set, 

AM 


&L 


/ \ 


/ 1 o \ 

1 — 1 
o 


z o 


i 



1 

II 

(N 

** 


\l 0 i 

1 

l 0 z ni 

\ ^ 


\ °l 


(60) 


We have seen by Eq. (6) that the quadratic form constructed from 
represents the power. From we obtain 

s( fe> * ?( J f^ + 8^- (l v+z o J l 2 + l V_Z 0 I | 2 ) 

(61) 


which is the sum of the powers in the forward and backward waves . 
In this form the conservation law is somewhat unfamiliar. But 
the net power s(Kjl) can be written also as the difference between 
the powers in the two waves , and we can form 
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I 


I 


K + 



(K n +K.) 

/V*X A^\Z 



K 







( 62 ) 


Then s(K+) and s (K^) are invariant also, and signify that power 
is conserved in each wave separately. 


We observe further that the metric 



0 


\ 


L 


/ 


( 63 ) 


where L and C are respectively the distributed inductance and 
capacitance per unit length, is also appropriate. Then 


= 


1 

2 


C V 


I L 


( 64 ) 


is invariant, and is recognized as the total electric and magnetic 
stored energy per unit length on the line. 

There may then be any number of invariants of physical inter- 
est. Only p can be independent, but these may not give explicitly 
the information we seek. We chose the metrics Hermitian, after 
all, in order to obtain an explicit physical interpretation. An 
examination of the space of Hermitian metrics with a view to 
identifying invariants of physical interest can lead also to a re- 
cognition of relationships that may not otherwise be obvious . The 
same relation will obtain among the s (K) as among the K. If 

A>v-\ MW 
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( 65 ) 


a.K. = 0 

AM 


then 


x 


I Y' a.K. )x = Y'' a.x^K.x a.s(K-j) 

I i—J iwilnw Z— / 1M* /•ni'WV Z — i 1 M*!' 

\i / i i 


= 0 


In this case. 


K , = K o * — (K + +K_) 

0^3 6l) AM 2 (jO {SjiX VA 


Then 


s ( K^s ) « [s (K+)+s (K_) ] 


AM AM 


( 66 ) 


(67) 


( 68 ) 


2 8 

The stored energy per unit length is thus times the sum of the 
powers in the forward and backward waves. 

As another example, a lumped shunt susceptance jb is re- 
presented on a voltage-current basis by the operator 


M 

/iaa 




(69) 


This operator has appropriate metrics of the form 


K 

/W\ 




(70) 


One choice of the arbitrary constants, f or , again gives the 

metric of Eq . (5), indicating that power is^conserved. Another 

choice gives 
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1 


0 





( 71 ) 


which implies that 

s (K„) = W* 

/WlZ 


(72) 


is conserved. The constancy of the voltage across the network 
appears here in a somewhat weakened form. 

The K of Eq. (70) can be chosen nonsingular, but K 2 is sin- 
gular (as^are K+ and K- in Eq . (60)). The space of appropriate 

Hermetian metrics wilT include singular metrics, even if we have 
gone to the trouble to choose the metrics that span this space 
nonsingular, and these may be of interest. The nonsingularity of 
the Khi does not ensure the nonsingularity of linear combinations. 
Conversely, a linear combination of singular K^i can be nonsin- 
gular (as in Eq. (67)). In attempting to interpret our results 
physically, then, it is not important that we choose Hermitian 
metrics by Eq . (59) to be nonsingular, only that we choose a 

linearly independent set. Later, however, it will be significant 
to know that a nonsingular set exists. 

The assumption that Hermitian metrics lead to quadratic 
forms that represent real physical invariants can also be mis- 
leading. For example, the operator for a lumped shunt resistance 
r , 



has the same repeated eigenvalue X = 1 as the operator of Eq. (69) 
for a lumped shunt susceptance jb. It will, therefore, also have 
appropriate Hermitian metrics, even though real power is not con- 
served. From one of these. 


K = 

/WV 




(74) 


with x as in Eq. (2) 

AAA 


we obtain 
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( 75 ) 


s (K) = i Im (VI*) 

rtV\ Z 

The imaginary part of the complex power — a real quantity, to be 
sure — is thus conserved. This example (see also Example 2 of 
Section VI) suggests that the special significance attributed to 
Hermetian metrics should not be pushed too far. 

While such simple systems are adequate vehicles for these 
several comments, examples of more interesting systems obeying 
conservation laws are given by Pease (refs. 1,2, and 4). 


Testing the Operator 


The conditions given for the existence of appropriate, non- 
singular K are conditions on the eigenvalues and structure of the 
system operator. It may be pointed out, however, that it is not 
necessary to obtain the eigenvalues, eigenvectors, and generalized 
eigenvectors directly, in order to determine whether these condi- 
tions are satisfied. 

The eigenvalues $ . of R are the roots of the characteristic 

, . 1 /W\ 

equation, 

| R- 3 1 1 = a B n + a , B n_1 + . . .+ a, 3 + a n = 0 (76) 

where I is the identity matrix. It is well known that the condi- 
tion for the roots of an algebraic equation to be real or occur 

in conjugate pairs is that the coefficients a^ be real, and so it 

is not necessary to solve Eq. (76) to obtain this information. In 
the case of the lumped-system operator, the eigenvalues Aj_ of M 
will satisfy Eq . (37) or Eq. (38) if and only if M and (M,t) - 1 ^ 

have the same eigenvalues , and thus the same chara / cteristic equa- 
tions, and it is only necessary to write these out and compare. 

This provides a preliminary check on a necessary existence condi- 
tion, relatively simple even for a matrix of large order. 

The suitability of the eigenvalues and structure of E. or M 
for the existence of nonsingular solutions to Eq. (8) or Eq. (33) 
can be determined by putting the A -matrices JXi - Al^ and R, 
or - AI and M ^ - AI, in their Smith normal forms. This is 
accomplished for a matrix defined over a field F, by a seq- 
uence of elementary A -matrix transformations defined by the follow- 
ing operations on A - AI: 

c /WA AM 
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1. The interchange of two rows (columns) . 

2. The multiplication of all elements in a row (column) by 
a nonzero scalar in F. 

3. The addition, to all elements of a row (column), of the 
corresponding elements of another row (column) multi- 
plied by a polynomial p(A) over F. 

By a systematic application of these transformations, A - AI can 
be put in its Smith normal form, /XM 



where each f^(A) is a polynomial in A, with leading coefficient 
1, which divides fi+i(A); where ±±±2 * * * fp / P = l,2,...,r, is 
the greatest common divisor of all p-square minors of A - A£; 
and where r is the rank of A. It can be shown that two matrices 
A, and can be related by a^similarity transformation if and only 
if their characteristic matrices A - AI and B - AI have the same 
Smith normal form (ref. 5) . ' v * 


A similar procedure, in which the complications of operating 
on polynomials are traded against the need to examine matrices of 
higher order, together with other procedures for obtaining infor- 
mation about the operator, will be described in Section III under 
"Testing the Kronecker Sum." The degree of f r (A), which can be 
shown to be the minimum polynomial of A (see Section IV) , gives 
some measure of the degeneracy in A, but not as much information; 
therefore, we will not pursue this^point. 


III. DEVELOPMENT BY KRONECKER-SUM FORMALISM 


In the previous section, we discussed conditions for an opera- 
tor to have appropriate nonsingular metrics, a formal procedure 
for determining the metrics, and the degeneracy of solutions. The 
metrics could also be obtained by writing out Eq. (8) or Eq. (36) 
in detail, and solving the resulting set of linear, homogeneous 
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equations. Short of solving these equations, we could determine 
the rank of the coefficient matrix and thereby the number of 
linearly independent solutions we can expect. In general, this 
procedure, applied to a set of homogeneous equations, would not 
tell us whether matrices formed in a prescribed manner from the 
vector solutions could be chosen nonsingular. Moreover, this 
would seem to be only an ad-hoc procedure — not necessarily simp- 
ler; and not to provide a formal theory of the conservation laws 
based on properties of the system operator. 

Writing out either of these equations, however, yields a co- 
efficient matrix that can be identified as a Kronecker sum (which 
incidentally provides an efficient means for obtaining this matrix) . 
The Kronecker sum has useful properties that will enable us to 
obtain expressions for the number of linearly independent metrics , 
in terms of the eigenvalues and structure of the system operator, 
by this approach also; and more significantly, to express the 
conditions on the operator for the existence of nonsingular solu- 
tions , and the degree and type of degeneracy in the operator , in 
terms of properties of a subspace of the space spanned by the 
eigenvectors and generalized eigenvectors of the coefficient 
matrix . 

It is the number of repeated zero eigenvalues and the number 
of linearly independent eigenvectors with zero eigenvalue, to- 
gether with certain upper bounds on these numbers , that must be 
determined in order to obtain the information we require. The 
upper bounds can be determined as the number of repeated zero 
eigenvalues and the number of linearly independent eigenvectors 
with zero eigenvalue of a related Kronecker sum in the operator. 
Matrices of the order n^ x n2 must be examined, but the subspace 
spanned by the eigenvectors and generalized eigenvectors with 
zero eigenvalue is the simplest to investigate. Moreover, the 
Kronecker sums can be compared under a simpler transformation 
than that to the Smith normal forms. In any event, a formal equi- 
valent to the theory developed in Section II is obtained by this 
approach, together with procedures for testing an operator that 
can be useful in particular problems and can, in fact, give some- 
what more information than the tests described in Section II under 
"Testing the Operator." The formal equivalent to the procedure 
for obtaining the metrics may be considered to be provided by the 
procedure to be described in Section VI. 


Kronecker Product 


The theory of the Kronecker or direct matrix product, and of 
the Kronecker sum, is a standard topic in matrix theory (refs. 3, 
6, and 7). The direct product is defined and some of its pro- 
perties are summarized below. 
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The direct product of two nxn matrices .A = (a-^j) and B = (b-^j), 
denoted by Axb, is defined as the n^xn2 matrix AM 

/VW AM 


Axb 

AM /W\ 


/ 


a ll~ 


a 01 B 

21am 


l 


nl 


B 


a 12^ ~ a lnS« 


a 22w 


\ 


- a 


B 


nn r 


■/ 


( 77 ) 


Here Axb is written as a partitioned matrix, with nxn blocks each 
consisting of the matrix B^ multiplied by an element a^j of A. For 
our purposes we restrict attention to square matrices .A and B of 
the same order, although the theory applies more generally. 


The elements of Axb can be written 

rtM AM 


(Axb) 

AM A M 


ij ,kh 


a 


ik b jh 


( 78 ) 


where double indices are used to indicate the ordering of the rows 
and columns of the direct product, each double index running 
through the values 11 , 12 , . . . , In , 21 , . . . , 2n, . . . ,nl , . . . ,nn, in that 
order . 

The direct product is distributive. 


(A+B) XC=AXC+BXC 

rwi /IM Ma AAA AM 

A x (B+C) =AXB+AXC 

^ JW\ AAA AAA AM M ^ 


and associative, 


A x (BxC) = (AXB) x C 

AAA AM AM AAA AM MA 


but not commutative , 


A x B B x A 

AAA AAA AM AAA. 

It also satisfies a matrix-product relation. 


(AxB) (CXD) = (AC) x (BD) 

AAA AM AM AAA AM AM AM 


( 79 ) 


( 80 ) 


( 81 ) 


( 82 ) 
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Let x = (x{ ) be an eigenvector of ^ with eigenvalue 

y(s)* , 'L (y|s) an eigenvector of B^ with eigenvalue y s . 

, u( r ' s ) as the n^-dimensional vector' having components 


A r , and 
Define 


(r , s) 
u . . ' 




(83) 


where the double index ij again represents a single index running 
through values 11 to nn in the same order as before. Then u( r » ,s ) 
is an eigenvector of A^xB^ with eigenvalue A r y g . rt * A 

If A and are semisimple, Eq. (83) gives n^ linearly in- 
dependent eigenvectors, and A*B is semisimple also. 

IW\ 

The eigenvalues of A*b, whether or not A and B are semisimple, 
multiplicative combinations A r y s , with all multipli- 
cities of the A r and the y s included. If A^ and B^ are nonsingular, 
their eigenvalues will be nonzero, and Axb will accordingly be 
nonsingular also. '*** 


The Kronecker Sum 


The Kronecker sum of two nxn matrices A and B, which we de- 
note by A (A, B) , is defined as the sum of two' direct products, of 

/u* w 

the form 


A (A,B) = A x I + I x B (84) 

t VAA f\AA KAA r\f*\ A* N AAA AAA 

where I is the nxn identity matrix. Like the direct product, this 
is not a commutative relation. 


A (A,B) ^ A (B , A) (85) 

A/VN AAA AAA AAA A*A 

If x( r ) is an eigenvector of A^ with eigenvalue Aj- , and y( s ) 
is an eigenvector of B with eigenvalue ys/ then u^ r ' s ^, with lw com- 
ponents defined by EqT - (83), is an eigenvector oT A (A ,B) with 

. — . AKA AAA XAA 

eigenvalue A r +y s . 

Aqain if A and B are semisimple, A(A,B) is semisimple also. 

nyv* AAA. At \ AAA AAA ftAA , 

In any case the eigenvalues of A (A ,BJ are the n z additive combina- 
tions A r +y s , with all multiplicities' of the A r and the y s included. 


Degeneracy in A(A,B) 

AAA AAA AAr\ 

From the eigenvalues A r of A and y s of B^, it is easy to deter- 
mine the number of repeated eigenvalues of A (A,B) equal to A r +y s , 

iw\ rwt ywv 
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for given values of r and s. There will be the same number of 
linearly independent eigenvectors and generalized eigenvectors of 
A^( A,B^ in chains with this eigenvalue, and these will span an 
eigensubspace . We shall have need to know the number of linearly 
independent eigenvectors in such an eigensubspace. In general, 
if A^ and are not semisimple, the derivation is involved, although 
the expressions are simple and will be given without proof. 

First, however, let us determine the dimensionality of a 
given eigensubspace. The indices r and s are assumed to run over 
values indexing the linearly independent eigenvectors, or separ- 
ate chains, of A^ and respectively. For each value of r there 
is a chain of A of length i r , and for each value of s a chain of 
B of length is? Then, each X r is of multiplicity £ r , and each 
ys is of multiplicity is • 


There can be more than one value of r for which the A r are 
equal, and more than one value of s for which the ys are equal. 
The number m r of repeated eigenvalues of A_ equal to A r , including 
multiplicities, is the same for each of these values of r. If we 
choose a representative value r from this set of r, we can denote 
the number of repeated eigenvalues of A^ equal to A r by m r , for 
each r in this set. Similarly, we denote the number of repeated 
eigenvalues of B equal to y~ by m^ - . Then we can write 

AAA ° 


E_ 


i 


r 


m— 

r 



s~s 


( 86 ) 


where by r~r we denote values of r such that A r =Ajr, and by s~s" 
values of s such that y s =y^r. 

By the property indicated in Section III under "The Kronecker 
Sum," the eigenvalue A r +y s of A,B) , for given values of r and s, 
will be repeated i r i s times. The number of repeated eigenvalues 
A r +y s , for all r~r and all s— s', is then 

r~r r~r s~iT 

s~s 


But there is yet another way in 
eigenvalues. For distinct values of 


nwn— (87) 


which A,(A r B) can have repeated 
r, r^r 2 / and of s, s"^^^ , 
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such that A— ^ A— and y— ^ y- , it can be that A— + y— = 
r l r 2 S 1 S 2 = *1 S 1 

A— + y— . Let us choose a representative pair (r,s) from each 
r 2 S 2 

set of ordered pairs (r,s") for which the A— + y— are equal. We 

denote the pairs in such a set by (r , s ) — (r , s) . Then the number of 
repeated eigenvalues A r +y s =Ap+ys, counting all r ~T and s—S - , for 
all (r , s ) — (r , s ) , is 


E . 

(r,s)- (r,s) 



r~r 


__E. 

(r,s)~(r,s) 


m-m— 
r s 


( 88 ) 


s~s 


Equation (88) gives also the number of linearly independent 
eigenvectors and generalized eigenvectors with eigenvalue AJ+yg. 

( IT S ) 

Among the eigenvectors will be the ii ' defined by Eq. (83) , 
for all r~r, s'— s', (r— s’) — (f ,s) . But unless either A^ or is semi- 

simple, there will be more. The number of linearly independent 
eigenvectors with eigenvalue A= + y= will be given by 

3T S 


£*rs 


E_E 

(r,s)-(r,s) r— r 
s— s’ 


min U r ,£ g ) 


(89) 


where by min (m,n) we denote the smaller of the integers m and n. 

Moreover, p== will be within the limits 
J. s 



min 


(r, s)~(r ,s) 




p=- < 
r rs 


E 


m-m— 
r s 


(r , s ) — (r , s ) 


(90) 


and will attain t|}e_upper limit if and only if & r = Z s = 1 for all 
r— r, s— s', (r , s") — (r , s ) . A necessary, but not sufficient, condition 
for pj^ to attain the lower limit is that there be only one r— r or 

only one s— s’ (i.e., SL— = m— or l— = m— ) for each (r,s") — (r,s) 

We shall be particularly interested in the special case 
where m^ = mg- for all (r,s") — (r,s) . The limits then become 
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(91) 


z 


m— < p== < 
r — r rs — 


£ 



where the summations are over all T for which there exist s such 

that (r ,s)~(r ,s) . In this case the lower limit_is attained if 
and only if £— = m— or £— = m— for each (r , s - ) ~ (r , s ) . 

3T 1C S S 

A significant feature of these results is that A(A,B^ can 
exhibit considerable degeneracy even when there is no degeneracy 
in A or B^. This will account for the fact that a nondegenerate 
operator can have a number of appropriate metrics . 


The Equation AX + XB T = C 

AMA« om mi am 


One application of the Kronecker sum is in the solution of 
the equation 


AX + XB T = C (92) 

«W> (W\ AM *AA (VM 


for = (xij), X an nxn matrix. Here A = (aj_j), B, = (b^j), and 

= (c-Lj) are given nxn matrices, and = (b j ^ ) is the transpose 
o'f B. By Eq. (78), and with the use of the 6 function. 


6 . . = 
1 3 


1 , i = j 
0 ,i ¥■ j 


Equation (92) can be written in detail as 


f a . , x, . tZ—zx^b 


E^-E,., 


k ik « ~“ ih “ jh kh 


S. , +6 . , b ., ) x, , 
ik jh ik jh kh 


-E 


(93) 


kh 


(Axl + IxB) . . , , x, , = c. . 

w w m l] ,kh kh 13 


If we define n -dimensional column vectors x" and c" having com- 
ponents x^j and c^j respectively, where in each case the double 
index is treated as a single index running over n^ values accord- 
ing to the ordering convention adopted at the beginning of this 
section, then Eqs . (93) are a rendering in detail of the single 
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matrix equation 


A (A, B) x" = c' (94) 

AAA <vxa vvv * 

We shall be interested, in this section, in the nontrivial 
solutions of Eq. (94) when c^' is the null vector. A necessary 
and sufficient condition for 


A (A, B) x" = 0 (95) 

AAA AJVV IWV (WV OAA 

to have a nontrivial solution is that the determinant of A^A,]^} 
vanish, or equivalently, that A^(A,B) have a zero eigenvalue, which 
may be repeated. If A^A^, B) has a zero eigenvalue repeated p m 
times, the rank r of can be anywhere in the range n 2 -p m <r< 

n 2 - 1. There will then be n 2 - r, or between 1 and p m , linearly 
independent solutions to Eq. (95) . The solutions xj' can also be 
thought of as vectors in the null space of A^A,!^) , the space 
spanned by the linearly independent eigenvectors of JY (A f> B^ with 
zero eigenvalue. This is a subspace of the zero - eigenvalue eigen- 
subspace of A (A,B) , the space spanned by the Pm linearly indepen- 
dent eigenvectors and generalized eigenvectors with zero eigen- 
value . 


These remarks apply in general to linear homogeneous equa- 
tions; but since A (A,B) is a Kronecker sum, we can be more parti- 
cular. The number ""of linearly independent engenvectors of A^A^/B^ 
with zero eigenvalue is given by Eq. (89) for + |ig = 0 , and is 
in the range indicated by Eq. (90) . 

Furthermore, in our application A^ and B^will be forms of the 
same operator. This will enable us to develop relations between 
Kronecker sums in different forms of the operator that are equi- 
valent to relations between forms of the operator described in 
Section II. 

Equation (92) has been studied by other approaches (refs. 8 
and 9) , and this is not the only route by which we could arrive 
at the results to be obtained in this section. But we shall be 
interested more generally in Eq. (92) for C = j(dX/dz), and we 
shall find that the operator in the form A7 a,B) can be identified 
with the conservation laws m a quite fundamental sense. Express- 
ing results in this section in terms of properties of the zero- 
eigenvalue eigensubspace of A (A,B) , then, exhibits them as one 
facet of the theory. 
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Structure of the Zero-Eigenvalue 

1* T 

Eigensubspace of A (R ,-R ) 

A/l* rJ* 

By Eqs . (92) and (94), Eg. (8) can be written in the form 


A(R + -R T )k" = 0 (96) 

/W <vm rW\ A/W aAA 

where k" is an n2-dimensional column vector formed from the ele- 
ments < o*f ^ in accordance with the convention adopted. The co- 
efficient matrix is of the form 


t T t T 

A (R ,-R ) = R xl-IxR= 

f\A a tvv\ r\nn rvn /iw iwi aaa 



(97) 

where the entries shown are nxn blocks . The nontrivial solutions 
to Eq. (96) are the eigenvectors of ^(R+,-R^) with zero eigenvalue. 

If we apply the results given in Section III (under "De- 
generacy in A(A,B)") in an analysis of A(R^,-R T ) under various 
conditions on the eigenvalues and structure of R^, we find that we 
can express the results of this analysis in terms of four signifi- 
cant parameters. *By (r , sj~ (r , s ) , we refer to all ordered pairs 
(r,¥) such that g^ = gg-, where T and S' index distinct eigenvalues 
of R^. The expressions in Section III (under "Degeneracy in j\^(A,B)") 
then apply with reference to the zero eigenvalue of A(R^,-RT). 
Otherwise the notation is as explained there. We note that for 
the eigenvalues of & to be all real or occur in conjugate pairs, 
it is required that the gr be all real or occur in conjugate pairs, 

and in addition that m— = m— for all (r,if)~ (r,s) . 

1C s 

The first parameter is the number of repeated zero eigen- 
values of A(R+,-R t ), which from Eq. (88) is given by 

✓Wl <W AM 
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(98) 


^m 


E 


m-m— 
r s 


(r,s)~ (r,s) 


where the summation is over all (r,s) such that 8— = f3g , and m— 
and m<r are the number of repeated eigenvalues g r and Bg of R, in- 
cluding multiplicities, equal to 8f and Bg- respectively. 

The number of linearly independent eigenvectors of A(R^,-R T ) 
with zero eigenvalue (the number of linearly independent metrics' 
appropriate to R) from Eg. (89) is given by 


P 


E_„ 

(r ,5-)~(r ,s) 



min 


r~r 


(H ,1 ) 

r ' s 


s~s 


(99) 


where we drop the subscript on p, since we will be concerned only 
with the zero eigenvalue. In Eq. (99), £ r and £ s are the lengths of 
the chains of R^ with eigenvalues B r and B s respectively. 

That mp and rrig in Eq. (98) , and £ r and £ s in Eq. (99) , can 
be interpreted as indicated, although strictly they refer to the 
eigenvalues and chains of R+ and -rT, is easily shown. 

am rvw 

The other two parameters are summations over all r, indepen- 
dent of the character of the eigenvalues of R. The first is de- 
fined by 


p ' 

^m 



( 100 ) 


Let the values of £ r for r~r be denoted by £ r , r~r, and be 
arranged such that r i 


£>£ , i— 1 , 2 ,... 

i l+l 

The fourth parameter is then defined by 

P' =^^(2i"l) * r< (101) 

r i 1 

r~r 
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The following relations can then be shown to hold: 


min < p<p m <p^ 


(r , s) ~(r , s) 


p<p' (102) 

n < P ' < P^ 


where 

A. p>0 if p m >0 , for then the set (r,s")Mr,s) is not empty. 
From another point of view, if p m >0 , A(R+,-R T ) is singu- 

AM >VA AM 

lar, in which case Eq. (96) has at least one nontrivial 
solution . 

B. p>n if the eigenvalues of R are real or occur in complex- 
conjugate pairs. In this case,p=n if and only if a single 
chain is associated with each distinct real eigenvalue 

of R and with at least one of each distinct conjugate 
pair"'. If chains of equal length are associated with dis- 
tinct conjugate pairs of eigenvalues, this condition re- 
quires that a single chain be associated with each dis- 
tinct eigenvalue of ,R. 

C. Pm = Pm if and only if the eigenvalues of R are real or 

occur in conjugate pairs. w 

D. p=p ’ if and only if the eigenvalues of R are real or 
occur in conjugate pairs, and chains of equal length are 
associated with conjugate pairs of eigenvalues. 

E. P ,= Pm if and only if R^ is semisimple. 

F. p'=n if and only if a single chain is associated with 
each distinct eigenvalue of R. 

Derivable from D and E, and from E and F, respectively, 

G. P=Pm if and only if the eigenvalues of .R are real or 
occur in conjugate pairs, and II is semisimple. 

H. p^=n if and only if the n eigenvalues of R^ are distinct. 

It follows from the above results that, if the eigenvalues 
of R are real or occur in conjugate pairs, 

n<p<p'<p m =p i ; (103) 
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The condition on the eigenvalues of R. , therefore, is sufficient 
to ensure that there will be at least n linearly independent solu- 
tions to Eq. (96) . There will also then be p' repeated zero eigen- 
values ^(R+,-R t ). But if R, and thus ^(R+,-gT) , is not semisimple, 
it is interesting to observe that A(R+,-RT) will have not only 
fewer than n 2 linearly independent^eigenTectors , but in particular 
fewer than p^ associated with the zero eigenvalue. It is charac- 
teristic of the class of results we have obtained that a property 
that depends on the eigenvalues and structure of a matrix depends 
only on the structure of the zero-eigenvalue eigensubspace of a 
properly formulated Kronecker sum. This will be illustrated more 
adequately as we proceed, after we introduce a slightly different 
Kronecker sum than the one we have been considering. We shall 
then find that the similarity requirement on R. and R ' for the 
existence of appropriate nonsingular metrics, which^is equivalent 
to conditions on the eigenvalues and structure of R, can be in 
like manner relaxed. w 


Reformulation of the Theory 

We have obtained explicit expressions for the number of lin- 
early independent metrics appropriate to a system operator JR, under 
various conditions on the operator, by an approach based on the 
Kronecker theory. We have expressed the results in terms of the 
parameters p, p m ,p', and p^, which are well defined in terms of 
the eigenvalues and structure of R. We have also established rela- 
tions among these parameters that^are equivalent to the conditions 
on the eigenvalues and structure of R, including the conditions 
under which the metrics can be chosen nonsingular. These results 
indicate a dependence of such interior properties of R., and the 
character as well as the number of its appropriate metrics, on 
conditions on p. Pm, p', and p^, . We have in effect reformulated 
our theoretical results in terms of conditions on these parameters , 
but the parameters are defined in terms of the interior properties 
of R. The unprimed parameters, however, are also well defined in 
terms of the structure of the zero-eigenvalue eigensubspace of 
A(R+,-R t ). We would like a comparable interpretation of p' and 
p^^in terms of simpler matrix properties. 

In attempting to characterize R and its appropriate metrics 
in terms of properties of A (R+ , -r t ) , n ~we encounter the difficulty 
that, while the parameter s < '*"p A *and~p m are descriptive of simple pro- 
perties of A (R+,-R^) , p' and p^ are values assumed by these para- 
meters under 1 conditions on R. We are able to bypass this diffi- 
culty, however, for we can show that p' and p^ are also descrip- 
tive of simple properties of a different Kronecker sum formed 
from _R. In fact, exactly as p and p m are defined in terms of the 
structure of the zero-eigenvalue eigensubspace of A^R^ , -R^) , p' 
and p^ are defined in terms of the structure of the zero-eigenvalue 
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eigensubspace of A(R,-R). 

/W*. HMV IW| 

Thus A (R, -R) will have 

am AM ran 


P ' 

^m 



( 104 ) 


repeated zero eigenvalues, where the summation is over all values 
of T, regardless of the character of the eigenvalues of R. The 
number of independent eigenvectors of A(R,-R) with zero eigenvalue 

• _ I# IWl «* 3 


P' 


= S X ( 2 i-i ) 

r i_ 
r~r 



( 105 ) 


where again the summation is over all r. The upper bounds on the 
dimensionalities of the zero-eigenvalue eigensubspace and the null 
space of A(R^,-rT), then, are respectively the dimensionalities 
of the zero-eigenvalue eigensubspace and the null space of A(R,-R). 

AM AM AAA 

We are thus able to restate the theory in the following 

terms : 


The number of linearly independent metrics appropriate 
to a system operator R. is given by the dimensionality 
of the null space of A(R+,-R T ), which is not greater 
than the dimens ionalilty^^nT of the null space of 
A (R, -R) . If the zero-eigenvalue eigensubspaces of 
X(R+-rT) and A.(R,-R), which contain the null spaces, 
are' of equal dimensionality, there will be at least 
n independent metrics. These can be chosen nonsin- 
gular if and only if the null spaces are also of equal 
dimensionality . 

By C and D, under "Structure of the Zero-Eigenvalue Eigen- 
subspace of A. (R+,-rT) » in Section III, only if the zero-eigenvalue 
eigensubspaces are^of equal dimensionality can the null spaces be 
of equal dimensionality. The former condition is necessary and 
the latter is necessary and sufficient to ensure that the metrics 
can be chosen nonsingular. 

We have also the following corollaries: 
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If the zero-eigenvalue eigensubspaces of A > (R+,-JR T ) and 
A (R,-R) are of equal dimensionality, and the null space 
of ^(R,-R) is of dimensionality n, then the null space 
of A ( S ? 1 ~' rT ) also be of dimensionality n, and the 

n linearly independent metrics can be chosen nonsingular. 

2. If the zero-eigenvalue eigensubspaces of A^(R^,-rT) and 
A(R, -R) are of equal dimensionality n, therT the* null 
spaces' will also be of dimensionality n, and the n lin- 
early independent metrics can be chosen nonsingular. 

Moreover, by C, E, and H: 

If the zero-eigenvalue eigensubspaces o f A > (R + ,-]R T ) and 
A(R,-R) are of equal dimensionality, and the null space 
of > ~A.(R,-R) is of this same dimensionality, then the null 
space of '"A (R + , -R T ) will also be of this dimensionality, 
and the metrics**can be chosen nonsingular. 

The condition given in Section II for the existence of non- 
singular metrics appropriate to R was that II and R+ be similar 
matrices , which required that the~eigenvaluesT and structure of R^ 
satisfy certain conditions . We have replaced these by the condi- 
tion that the null spaces of A(R^,-R t ) and A(R,-R) be of the same 
dimensionality (p=p 1 ) . This condition can be expressed also in 
terms of the rank r=n 2 -p of A(R+,-R T ) and the rank r'=n 2 -p' of 
A (R , -R) , and is that w A *' w 

fl/w nnA am 

r = r' (106) 

or that A(Rt-R T ) and A(R,-R) have the same rank. 

AM AM AAA AAA AM 

It is well known that two matrics have the same rank if and 
only if they are equivalent ♦ (By definition, A and are equiva- 
lent if there exist nonsingular matrices P and Q such that A=PBQ.) 
Thus at the sacrifice of an increase m the order of the matrices, 
we have replaced the requirement that and R+ be similar (which 
requires that Q=P“1) with the weaker requirement that A(R+,-R T ) 
and A (R , -R) be equivalent, 

rsAA AAA 

We shall encounter, in the test for semisimplicity described 
below, another example of an interior, structural property of R 
that can be expressed as the equivalence of two Kronecker sums . 


Testing the Kronecker Sum 

These results provide a basis for testing an operator to 
obtain information about its appropriate metrics — and indeed about 
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the eigenvalues and structure of any matrix, apart from our inter- 
est in the conservation laws governing a system it might describe 
— by determining or comparing the dimensionalities of the null 
spaces and zero-eigenvalue eigensubspaces of appropriate Kronecker 
sums . 


Test I : test for suitability of eigenvalues of By C, 

in Section III, under "Structure of the Zero-Eigenvalue Eigen- 
subspace of A (R+ , -rT) , " the eigenvalues of R are real or occur 
in complex-conjugate pairs if and only if Pm=Pi{i. These parameters 
could be determined as in Test IV below, and compared. We indi- 
cate this approach only for consistency, as a formal alternative 
to the much simpler test described in Section II under "Testing 
the Operator." That preliminary test, which requires only writing 
out the characteristic polynomial of R, is clearly indicated in 
any test procedure. 

Test II ; determination of the number p of linearly indepen- 
dent metrics appropriate to R . - Since p is also the dimensional- 
ity of the null space of X ( Rt , — ) , it can be found from the rank 
r of A (R+ ,-RT) , for r=n.2-p . *** The rank r can be determined , among 

/MA /VV* . 

other ways, by a sequence of elementary transformations defined 
by the following operations: 

1. The interchange of two rows (columns) . 

2. The multiplication of all elements in a row (column) 
by a nonzero scalar. 

3. The addition, to all elements of a row (column) , of the 
corresponding elements of another row (column) multi- 
plied by a scalar. 

By a systematic application of these elementary transformations, 
any matrix can be put in its normal form . 



where Ij- is the rxr identity matrix (ref. 5). We obtain p, then, 
from the order of I r . 

Test III: test for suitability of eigenvalues and structure 

of - By D, the eigenvalues of R are real or occur in conjugate 
pairs, and chains of equal length^are associated with conjugate 
pairs of eigenvalues, if and only if p=p ' — that is, A(R^,-R^) and 
A(R^,-R) are equivalent matrices, having the same rankT ^The ranks 
cTf'^these matrices can be determined, and compared, by putting each 
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in its normal form. Unlike putting R-f3I and R^-$I in their Smith 
normal forms, this does not involve the manxpulation of poly- 
nomials and the determination of greatest common polynomial divi- 
sors. Furthermore, it is a property of equivalent matrices that 
one can be obtained from the other by elementary transformations, 
which permits a comparison without reducing each to its normal 
form. 

Test IV : test for se misimplicity of^.- By E, R is semi- 

simple if and only if p'=Pm- Since p^ is the number of repeated 
zero-eigenvalues of A(R, -jR), in principle p^ could be determined 
by finding the term of**lowest degree in the characteristic poly- 
nomial of A (R^, -R) for which the coefficient does not vanish. Thus, 
a measure of the number of repeated eigenvalues of .R, determined 
conventionally by finding the roots of the characteristic poly- 
nomial, can be obtained merely by inspection of the characteristic 
polynomial of A (R^, -R^) . But while not requiring the solution of 
an algebraic equation, writing out the characteristic polynomial 
(or even determining a sufficient number of the polynomial co- 
efficients to apply this test) involves considerable, if straight- 
forward, calculation for a matrix of high order. Compared with 
the tests we are describing on n 2 xn 2 matrices, it is relatively 
simple to do for an nxn matrix, and is therefore indicated above 
in the first of these tests . But rather than write out the char- 
acteristic polynomial for an n 2 xn 2 matrix, we can avail ourselves 
of a more expeditious procedure to test R^ for semisimplicity. 

Denote the number p' of independent eigenvectors of A(R,-R) 
with zero eigenvalue by p'(A). In all there are p^ linearTy^in^ 
dependent eigenvectors and generalized eigenvectors with zero 
eigenvalue, in chains x^'^', x ' 2 ) , . . .x vi .&i) such that 

AM 


A (R,-R)x (l/1 ^ = 

0 
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A (R,-R)x (l,2) = 

xU, 1 ) 

aaai a/v» Am va 

AM 

P* 

j. 

(i,£.-l) 

A (R, -R) x = 

AAA. AAA MA AM 

x 1 

MA 


From the first two of these equations, we see that 
[A (R , -R) ] 2 x* 1 ' 1 * = 0 

AM AAA Mn vTto An* 

[A (R, -R) ] 2 x (l ' 2) = A (R,-R)x (l,1) = 0 

am Am AM M A. AM AAA AM. AM 
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Thus the x (i,1) and x^' 2 ) are eigenvectors of [A(R,-R)] 2 with 

. -V* AM AM AM. AM 

zero eigenvalue. 

2 . 2 
There are then p 1 (A ) independent eigenvectors of [A^R^-R)] 

with zero eigenvalue , consisting of the eigenvectors and general- 
ized eigenvectors of rank 2, with zero eigenvalue, of A^R^-R). If 
the chains are not all of length 1, p 1 (K)<p 1 (A/) • Similarly, 
there are p'(A3) independent eigenvectors of [A(R,-R)]3 with zero 

am am AAfc AM , , 

eigenvalue, consisting of the eigenvectors and generalized eigen- 
vectors up to rank 3, with zero eigenvalue, of AJR,-R) ; and if 
the chains are not all of length less than 3, p' (A2)<p' (A 2 ) . 
Finally, if £ m if the length of the longest chain v 'associated with 
the zero eigenvalue of TV^R^R) , then the null space of [ AjR^ -R) ] 
and of all higher powers" of 2£(R/“R) , consists of the entire p^ - 
dimensional zero-eigenvalue ’eigensubspace of A(R,-R). We there- 

_r- , AAA /\AA AAA 

fore have 


n < p* (A) < p 1 (A ) <, 

AAA 


Z Z + 1 

,< p’ (A m ) = p’ (A m ) = p* (107) 
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or denoting the rank of [A(R,-R)]^ by r' (A*^) , 

AM AM AM AAA 

n - n > r' (A) > r' (A 2 ) > . . . > r ' (A m ) = r ' (A m ) = n - p' 
am am ma /w% m 

(108) 

since r' (A^) = n 2 - p ' (A^) . 

/IM AAA 

Since by Eq . (107) p' (A)=p^ requires that p' ( A ) =p ' (A^ ) , it 

follows that R is semisimple^if and only if A (R, -R^and (A(R,-R)]^ 
are equivalent, having the same rank. The ranks of these matrices 
can be compared as described in the previous test. 

The value of p^ could be found by determining the rank of 
the smallest power of A^R^ - **) having the same rank as the next 
higher power. ( Similarly", ^=n 2 -r(Au m ) could be found by deter- 
mining the rank r (A^ m ) of the smallest power of ^(Rj,-Rj) having 
the same rank as the next higher power; but, if & m =l, it does not 
follow that R is semisimple.) 

Test V : test for multiplicity of chains associated with 

distinct eigenvalues of & . - By F, there will be a single chain 
associated with each distinct eigenvalue of ^ if and only if p'=n. 
As described in Test II, p'=n 2 -r' can be determined from the rank 
r ' of A (R, -R) . 

Ann haa fwi 

Test VI : test for distinct eigenvalues of & . - By H, the 

eigenvalues of R are distinct if and only if Pm =n f and thus by 

*AA 
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Eqs. (107) and (108), if and only if [A(R,-R)] 2 is of rank n 2 -n. 

«/VI AM 

These tests can be applied independently, each assuming no 
special conditions on JR that must be determined by another test. 
But as we would expect from the theory on which these tests are 
based, results of some tests can in certain instances be combined 
to infer the results of others. If, for example, it is found by 
Test I that the eigenvalues of R are real or occur in conjugate 
pairs, then Eq . (103) is valid.*''’’ If it is further determined, 
either by Test V that p'=n, or by Test VI that p r [ t =n, it follows 
that p=p'=n. Then without having determined p directly, we know 
that there are n independent metrics appropriate to ^R, and with- 
out Test III that these can be chosen nonsingular. Or if by Tests 
I and IV it is determined that the eigenvalues of ^R are real or 
occur in conjugate pairs, and .R is semisimple, then by E and G, 
p=p'=p I J 1 >n; and again it is not necessary to apply Test III in 
order to conclude that the metrics can be chosen nonsingular. 


Lumped Systems 

The theory concerning the metrics appropriate to a lumped 
system can be reformulated along the same lines that we have dis- 
cussed in detail for the uniform distributed system. Using the 
technique described in Section III under "The Equation _AX+XB T =C, " 
we can write Eq . (36) , the condition for a metric K to be approp- 

riate to a lumped-system operator M, in the form 

OAA 

A (M + ,- (M T ) _1 ) k" = 0 (109) 

am rw\ av» rtAA nAA 

The relations among the properties of AJM^ , - (iM T ) ) and A(M,-M), 
the eigenvalues and structure of M, and the number and character 
of the nontrivial solutions to EqT (109), are then in complete 
parallel with those among the properties of yV^Rt,-!^) and ^(R^-R), 
the eigenvalues and structure of .R, and the number and character 
of the nontrivial solutions to Eq . (96) — where again p and p' 

are the dimensionalities of the null spaces, and p m and p^ of the 
zero-eigenvalue eigensubspaces, of A (Mt ,-(M T ) ~1) and A^M, -M) res- 
pectively. It is only necessary to" substitute everywhere, for the 
conditions given by Eqs. (12) and (13) on the eigenvalues of R, 
the conditions of Eqs. (37) and (38) on the eigenvalues of M.** 

AAA. 

Alternatively, we could first obtain an equivalent .R for the 
given M, as discussed under "Lumped Systems" in Section II, and 
reduce the problem to the form we discussed in detail. 

We have thus far been occupied with a parallel treatment of 
the operators R^ and M, or jM(z,zi) , for constant R., in terms of 
equivalent conditions on their eigenvalues and structure, translated 
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into identical conditions on certain Kronecker sums, and with 
reference to the same set of constant metrics. We shall extend 
this parallel through the following section. Later we shall 
exploit some dissimilar properties of JR and M.(z,zp), and intro- 
duce other metrics whose dependence on the operators does not 
exhibit such simple parallels. 


XV. THE CONVERSE PROBLEM 
Distributed Systems 

In the previous sections we addressed ourselves to the prob- 
lem of what can be said about the metrics appropriate to a system 
operator, for distributed systems of a restricted class, and those 
appropriate to a lumped-system operator. Some attention has also 
been given to the converse problem, for uniform distributed systems 
of how an operator is determined by its appropriate metrics (ref.l) 
A theoretical treatment of this problem is of interest in that it 
permits the synthesis of systems governed by specified conserva- 
tion laws, and more basically in that it points to the significance 
of the conservation laws in defining a system. A critical look 
at the published results will show that they are only partly satis- 
factory in providing such an alternative approach to system defini- 
tion. (See also Section VII.) 

Given any two independent, nonsingular, Hermitian, constant 
metrics Kj_ and K 2 , the range of operators R^ for which and £2 
are appropriate metrics can be determined. For Kp to be appropri- 
ate to R, by Eq. (8) it is necessary and sufficient that KpR^ be 
Hermetian — i.e., that 


K,R = H 

«wJL mw am 

where H is any Hermitian matrix. Thus 

R = k7 x h (110) 

am m*-L am 

gives the range of JR having K-^ as an appropriate metric. For £2 
to be appropriate also, it is necessary and sufficient that 

R + K„ = K„R 

where R^ is in the range given by Eq. (110) . On substitution of 
Eq. (110) , this relation can be put in the form 
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J + H = HJ 

ATf\ AW\ Alft 


where J is defined by 
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( 112 ) 


Equation (111) places a further 
tor R of the form of Eq. (110) , 
to Eq\ (111) , will have both K, 


restriction on H,. Thus any opera- 
where H is an Hermitian solution 

_ /W\ 

and K as appropriate metrics. 


Equation (110) defines a class of systems with operators 
having as an appropriate metric, such that s(K>l) is an invari- 
ant. Similarly, R = K^H, where H is any Hermetian matrix, de- 
fines a class of systems^for whicli s(I< 2 ) is invariant. The class 
of systems for which both s(Kq) and s(K2) are invariant can be 
defined as the intersection of these two classes. 


The extent of this intersection, of the range of .R , will de- 
pend on the choice of JKq and K 2 / for these may imply other indepen- 
dent metrics that must be appropriate to any R in this range. It 
is easy to verify that if Kg and K 2 are appropriate to R, then 
every metric of the set w *** 


K 

f*n m 


= K n J 

(Wnl W 


m-1 


m=l , 2 , . . . 


(113) 


is also appropriate. For by Eqs . (110) to (113), 

R f K = HK7 1 K 1 J 01 ”^ = HJ 1 "- 1 = J f HJ m ' 2 = ... = (J + ) m_1 H 


K R = K-, j m 1 K 1 1 H = K.j m 2 K. 1 K„K 1 1 H = K-, j m 2 R n 1 J + H 
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= K-, K-, 1 (J + ) m 1 H 

1 •*\J_ /+* A*** 
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Since every matrix satisfies its characteristic equation, the n*n 
matrices J m_ l, for m>n, must be linear combinations of the J m- 1, 
m=l, 2,./T,n. Thus at most K]_ , K 2 / • • • / K n (and at least JKp , K 2 T will 
be linearly independent. The more independent me trics ,,, ’ implied by 
a specified and K 2 , the more restricted will be the class of 
systems that can be designed that will be governed by the conserva- 
tion laws ds(K^)/dz = 0 and ds(K 2 )/dz = 0. 
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Specifying additional metrics, other than those implied, 
will further restrict the class of systems — ultimately to R=I, 
for which all metrics are appropriate. It is in this sense tfTat 
a system is said to be defined by its appropriate metrics.* 


Let us, however, look at these results in a little more de- 
tail. Let us first ascertain the minimum range of ^ for specified 
Ki and K 2 . With defined by Eq. (112) , and recalling the assump- 
tions on and , it is easily checked that 


H 

(W 


K 1 


H 

AM 


&2 


(114) 


are independent, nonsingular solutions to Eq . (111). But Eq.(lll) 

is of the form of Eq. (8) , and we saw in Section II that if an 
equation of this form has one nonsingular solution, it has at 
least n independent solutions, which can be chosen Hermitian .Thus , 
by Eq. (110) there will be at least n independent operators R 
having Ki and K 2 as appropriate metrics — regardless of the number 
of additional "independent metrics KiJj 11- -*-, m=3, 4,..., that are 
implied by K, and K 0 , and that restrict the range of R. 

The number of independent metrics KiJ^ ^ is inversely related 
to the degeneracy in J. In the notation of Section III, if differ- 
ent subscripts are used to index repeated eigenvalues of J when 
these are associated with independent eigenvectors or separate 
chains, and if £ r is the length of the chain with eigenvalue r\ T , 
the characteristic equation of J can be written 
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(n-n r 
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(115) 


and is of degree n • 
tion of lowest degree 
form 


However, this may not be the polynomial equa- 

that J satisfies. The latter will be of the 

«*> 


max 

|J (n-n-) 

r 



r~r 


0 


(116) 


where the n r are distinct, and max £ r , r~ r, is the length of the 
longest chain with eigenvalue rir = ^r* The polynomial on the left 


The properties of an operator on a wave basis having an approp- 
riate metric of canonical form for pairwise coupling have also 
been arrived at by an approach beginning with the specification of 
the metric (see ref. 10). 
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in Eq. (116) is known as the minimum polynomial of ^J, and is of 
degree not exceeding n. The less degeneracy in .J, The larger 
will be the degree of the minimum polynomial, and thus the greater 
the number of linearly independent metrics 

But applying the results of Sections II and III, we see that 
the less degeneracy in the fewer linearly independent solutions 
H there will be to Eq . (Ill), and the smaller will be the range 
o*f R. This explains the mechanism by which the range of R is re- 
stricted by the specification of metrics and J&2 that imply addi- 
tional metrics. The degeneracy in will be minimum when J is 
semisimple with distinct eigenvalues, or when only one chain is 
associated with each distinct eigenvalue. The minimum polynomial 
will then be the same as the characteristic polynomial, and there 
will be n independent metrics in the set But we have 

seen that under this condition, and only Then, will there be only 
n independent solutions to Eq. (Ill) , and thus a minimum n-dimen- 
sional range for R. 

rwi 

We also note that 


R, - J k 1 r k = 1,2, .. . (117) 

rm K f\M 


are operators in the range of I*. For 
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r t x k-l 


1\k-2 „ -1 T 


(J x ) ' K, = (J ) K, = (J ) K-K. K-, = (J') 
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= (J + ) k 2 K, J 
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= ... = K-, J 
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(J k 1 ) t K = (J + ) k-1 K = (J + ) k 2 K K 1 K 
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(J f ) k - 2 K,J 
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. . . = K 0 J 
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k-1 


and so K n and K„ are appropriate to each R. . 

rvM_L AM Z 

If the degeneracy in J is minimum, there will be n indepen- 
dent R k . But the range of R is then only n-dimensional , and so 
there**are no more independent operators in this range. Further- 
more, there will then be n independent metrics K =K,J m- l, m=l,2,...,n 
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which are all appropriate to each of the JR^, k=l, 2 ,...,n. 

We cannot in general further restrict the range of R. selec- 
tively by the specification of additional metrics, such as to de- 
fine each independent It in this range uniquely in terms of its 
appropriate metrics . Any metric will be appropriate to Ri=^. 

Since the degeneracy in £ is minimum, R2 = ^, will have exactly n 
independent, appropriate metrics; i.e.T only the JK m that we have 
already determined. Further degeneracy may or may not appear in 
successive powers of J. For example, if R2 = ^. i- s semisimple with 
distinct eigenvalues rir/ then Rj c =j| c- 1, k=3, 4,..., will be semi- 
simple with eigenvalues • '"’clearly the ri£-l, for any k> 2 , may 

or may not be distinct (e.g., if rii=l, ri 2 =- l' then r\^=r\2=l; but a 
situation of this kind is not representative) . If successive 
powers of J do not become more degenerate, the m=l,2,...,n, 

will be a maximal set of independent metrics appropriate to each 
of the independent operators R k , k=2,3,...,n. (It is easily 
verified that another independent operator can be included, making 
a total of n, by adding to any one of the JJk, k=2 ,3, . . . ,n.) 

If the degeneracy in J is not minimum, the range of R. will 
be greater than n-dimensiorial . There will still be q independent 
operators in this range of the form Rk = ;j] c_ '*' ' k=l, 2 ,...,q ( 2 <q<n), 
where q is the degree of the minimum polynomial of J^; and the q 
independent metrics JK m =Kp J* 11 - 1 , m=l,2,...,q, will be appropriate 
to each of these operators’. Since the degeneracy in JT is not 
minimum, the total number of independent metrics appropriate to 
R 2 =J will, however, be greater than n. But let Kf be a metric, 
"appTopriate to R 2 , that is independent of the K m , m=l, 2 ,...,q. 

Then 

J + K. = K.J (118) 

/VV) /vwl ** 1 ' m 

Postmultiplying both sides by J, and making use of Eq . (118) in 

the result, we obtain A 

(J 2 ) + K. = K.J 2 (119) 

«*r w«l <W> 

and, accordingly, is appropriate also to ^ 3 =^ . By extension, 
Kf, and thus all o'? the metrics appropriate to R 2 , are appropriate 
to each Rk, k=l,2,...,q. Again, further degeneracy may or may not 
appear in successive powers of J. (Consider, for example, a semi- 
simple J with repeated eigenvalues. The distinct eigenvalues 
raised To a given power will normally, but not necessarily, be dis- 
tinct.) If successive powers of <J do not become more degenerate, 
the independent operators Jlj., k=2T3,...,q, (and R 2 +;E) will have 
the same set of appropriate metrics. 
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We can easily construct examples also of nonsemisimple J, 
for which the degeneracy does not increase in powers of And 

if the degeneracy in J is not minimum, we can often find a set 
of operators in the range of R that are independent of the Rj. , 
k=l,2,...,q, and which also have a common set of appropriate 
metrics. In fact, with R substituted for Eq. (119) follows 
from Eq. (118) for any R in the range; and while the degree of 
degeneracy in such an R need not be the same as that in J, inde- 
pendent powers of R having the same degree of degeneracy will 
have the same appropriate metrics. 

We have indicated how the degeneracy in determines both 
the range of £ and the number of metrics implied by JKi and £ 2 , 
and gives some measure of the extent to which this range can be 
subdivided by the specification of additional metrics. We have 
found that specifying a set of nonsingular metrics, and stipula- 
ting that it be maximal, still in general determines R only with- 
in a range. 

Furthermore, this approach to system definition, at least in 
the form given, lacks generality. The procedure does not apply 
to a specified set of metrics that cannot be chosen nonsingular 
in linear combination, but which also define a range of operators. 
And there is a range of uniform systems that do not obey any con- 
servation law. Specifying the constant appropriate to a non- 
uniform-system operator would in general be even less useful for 
discriminating among systems. (It does not follow from the exis- 
tence of one constant, nonsingular solution to Eq. (8), or Eq.(36) 
for M a matricant, that there need be more, if R is not constant 

<W\ fVY\ 

— but we shall see in Section V that the specification of z-de- 
pendent solutions to the latter equation will lead to a different 
interpretation of the conservation laws governing a range of sys- 
tems determined by a variation of this procedure.) The results 
of this section suggest that the constant K might play a role in 
a more general, exact scheme for system definition; but clearly 
something more must be incorporated in this scheme. What this 
must be, and its physical interpretation, will become clearer as 
we proceed, and the question will be resolved finally in Section 
VII. 


Lumped Systems 


We next ask how the range of lumped-system operators having 
specified appropriate metrics can be determined. 

A set of constant K defines a range of constant R, and there- 
by a range of matricant’s as solutions to Eqs . (39) and 1 (41) for R 
in this range. A range of lumped-system operators M that are K- ,w 
Unitary with respect to every K in this set is giverT by the range 
of matricants evaluated at = 0 and any z = L. 
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The metrics of Eqs . (5) and (71), for example, define con- 

stant R in the range 

AM 


R 




( 120 ) 


where a-^ and a 2 are real constants. The range of matricants for 
R in this range is given by 

MO 


M (z , 0) 




( 121 ) 


Evaluated at z = L ^ 0 , with a± = 2irn/L and a 2 = b/L (n=any in- 
teger) , this is the operator 1M of Eq. (69) for a lumped shunt 
susceptance jb, and for which we found the specified metrics were 
appropriate. More generally M describes a lossless system at 
the edge of a pass band - reciprocal and symmetrical if 
e“ j a 1 L = + 1. 

It is interesting to note, as in this example, that certain 
metrics may define distributed systems that do not appear to be 
of physical interest, but which in turn give rise to lumped 
systems of interest. 

For a given R, M(L,0) depends on the choice of L. Further- 
more, K-Unitary M^carT be obtained as the matricant, evaluated 
at z =,W L, of JR outside this range, including R(z) , and can des- 
cribe a variety of physical systems. 

The significant point, however, is that for any choice of 
z = L 7 ^ 0 , the technique described gives the complete range of 
nonsingular M. for which the specified K are appropriate. For 
given any M Tn this range — which couTd describe a section of 
a nonuniform distributed system, a network of lumped elements, 
or an arbitrary transformation — £1 has an equivalent representa- 
tion as a section, of length L, of a uniform system described by 
an operator R for which these metrics are appropriate. This 
operator wilT therefore be in the range of constant R defined 
by the specified metrics, and M will be its matricant evaluated 


Even for a given choice of z = L, there will not be a one- 
to-one correspondence between iM and R^ in their respective ranges . 
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As indicated in Section II under "Lumped Systems," and as may 
be seen in the above example, there are still an infinite number 
of choices for R in its range from which we can obtain a given M. 

But the range of M, which is independent of the choice of 
L, is determined by the rang e of constant R defined by the speci- 
fied metrics. Lumped systems can therefore be treated by the 
techniques developed for uniform distributed systems. The rela- 
tion between operators for lumped and uniform distributed systems 
will be indicated more formally in Section VII, and it is parti- 
cularly appropriate that they can be so related. For if we 
further delimit the range of JR, in a manner to be discussed, we 
will in general further restrict the range of £1; but unless this 
is done by specifying additional constant JK, we shall not be 
distinguishing between operators in this range on the basis of 
what can be interpreted strictly as conservation laws. We shall 
find that this is a feature of uniform systems in particular. 


V. NONUNIFORM SYSTEMS 

In this section we generalize some of the previous results 
to nonuniform systems, and we introduce some metrics that we 
allow to vary with z in a manner appropriate to the conservation 
laws governing such systems. We consider two ways in which this 
can be done, each proceeding from a different interpretation of 
the conservation laws and providing different insights into the 
more general problem, but either providing a basis for analysis 
or synthesis. By ascertaining what is required of a metric for 
it to determine a conservation law, we justify the assumption 
that these must be constant for uniform systems as a special case. 
But the z-dependent metrics introduced in this section will later 
have application in characterizing uniform systems as well. 

The condition for a quadratic form s (K) , .K constant, to be 
invariant under the operation R, was given in Section II to be 
that K satisfy 

am 

R + K = KR (8) 

/vn AM AM AM 


This condition was derived without assuming R constant, and 
accordingly it applies to both uniform and nonuniform systems. 

But the subsequent conditions on R, or on A (R^ , -RT) and 
A (R,-R), for the existence of constant (nonsingular or singular) 
solutions to Eq . (8), the expressions for the number of these that 

will be linearly independent, and the procedures for obtaining 
constant solutions were arrived at on the assumption that is 
constant . 
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If R is not constant, the same procedures can be used 
formally^to obtain solutions to Eq. (8) . The number of indepen- 
dent solutions will depend as before on the conjugacy properties 
of the eigenvalues of R, the dimensionality of the eigensub- 
spaces , and their organization into chains of generalized eigen- 
vectors — any or all of which, however, may depend on position. 
But if they do not, then on a Jordan basis they will determine 
a number p of independent, constant solutions to Eq. (19) , 


,+ 

AM 


K' 


K ' R 1 


(19) 


where p has the same dependency as before on the eigenvalues and 
structure of — although the eigenvalues, eigenvectors, and 
generalized eigenvectors are allowed to vary with z subject to 
the above constraints. If the eigenvectors and generalized 
eigenvectors, and thus the modal matrix S^, are constant, the 
K' can be transformed back to the same number p of independent, 
constant metrics K on the original basis by Eq. (24) . In this 
case, all of the Theoretical results given for constant 11 apply 
(although, as we shall see, we must be somewhat more careful in 
interpreting the conservation laws) . The theory developed thus 
far, then, applies to a somewhat broader class of operators 
than assumed. If the modal matrix is not constant, however, 
there is no guarantee that any or all of the K' , in any set of 
independent situations to Eq. (19), will transform to constant K. 

AW 


There will still be the same number p of independent solu- 
tions to Eq . (8) . But if IC = .K(z) , then apart from any question 

of s (£) having the same interpretation for all z in terms of 
properties of the system, Eq. (8) no longer gives the condition 
on K for s (K) to be invariant. Equation (7) then becomes 

rw\ /w\ 


ds = d^_ 
dz dz 


Kx + x + 
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and the condition for s (K) to be invariant is more generally 

A*A 


R + K - KR = j 

ru* /w> n CLZ 


( 122 ) 


For constant K, Eq. 


(122) reduces to Eq . (8). 
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( 123 ) 


If s (K) 

<WV 


+ 

X 

nm 


is invariant, and K = K(z) , then 

rwv Avi 

(z)K(z)x(z) = x + (z, ) K (z, )x (z, ) 

yw\ am am X X am X 


where zj is any fixed value of z. By Eq. (42), the matricant 
M(z,zi) must then satisfy 

AAA 

(z , z, ) K (z)M (z , z, ) = K (z, ) (124) 

(V*> X aao am X VIA X 

Equations (122) and (124) are equivalent, and one can be obtained 
from the other with the use of Eqs . (39) and (41). 

Equation (122) is of the same form as Eq. (92) , with C = 
j(dK/dz) chosen as the derivative of the unknown JK instead of as 
a constant matrix. The derivation of Eq . (94) from Eq. (92) is 

not affected if is of this form, and accordingly Eq. (122) can 

be written 


A(R + ,-R T )k" = j (125) 

OM AAA AV> LIZ. 

This is a differential equation in the form of Eq. (1) , where in 
this case the operator A(R+,-R T ) is a Kronecker sum, and k" is an 
__ -dimensional state vector formed from the elements of £ in the 
manner described in Section III. We shall find in Section VII 
that k" is a state vector of the system in not only a formal, but 
a very real sense. If R satisfies the Lipschitz condition, 
A(R+,- rT ) will also? for 'its elements are simply elements of -R > , 
'con'juga'tes of elements of R, or sums of these. The same condition 
that ensures the existence^of a matricant, then, ensures the 
existence of n^ independent solutions to Eq. (125). The results 
of Section III give, for R constant, the number p of independent 
solutions that can be chosen constant. 

Thus it is only necessary that ^ satisfy the Lipschitz con- 
dition for there to exist n^ independent, invariant forms s (Jp . 

Let us look more carefully, for the general case where R is not 
restricted to be a constant operator, at the conditions for such 
forms to represent invariant properties of the system. A system 
is defined by the elements rij of its operator R, and the vari- 
ables chosen as components of the state vector jc. A consistent 
definition of what we mean by any property of the system requires 
that it depend in the same way on the rij and the components of 
x everywhere in the system — that is , that it does not depend 
explicitly on z. Thus it is reasonable to interpret an invariant 
form s (K) as a conservation law, although K is a function of z, 
if K carT be written as an explicit function' of the rij (z) only: 
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( 126 ) 


K = K (r . . (z) ) 

r** AAA -L J 

Constant K are a special case of Eq . (126), for both uniform 

and nonuniform systems. If the r-^j are constant, then conversely 
must be constant to be of this form. Thus we are justified in 
having restricted our attention to constant IC in considering the 
conservation laws governing a uniform system. But for a nonuni- 
form system, any solution K^(z) to Eq. (122) will meet this require- 
ment if, on substitution of z expressed in terms of any of the rj_j 
by inverting rij = rj_j (z) , K^rij) is single valued. 

All s (K) = f(xi,rij) derived in this manner from z-dependent 
solutions to' Eq . (122) are not likely to be identifiable as famil- 

iar invariants such as the power or the energy. But there is no 
guarantee either that the p independent invariants we have been 
considering for uniform systems will have simple physical inter- 
pretations, and Eq. (126) is the logical extension of constant 
to nonuniform systems. Furthermore, if the system is nonuniform, 
it is only by an examination of this larger class of metrics that 
we can be sure we have identified all the invariants of clear 
physical interest. (See example in the following section.) 

Although, for constant K^, an invariant form s (K) always re- 
presents an invariant property of the system, by the same token 
we must be careful in interpreting this property if the system is 
nonuniform. Suppose, for example, that the system includes a 
capacitance that varies in some manner C = Cog(z), g(0) = 1, and 
it is determined that (1/2 )Cq | V | 2 is invariant, where V is the volt- 
age. If the system were uniform, g(z) constant, we could conclude 
that the stored electric energy is conserved; but in general we 
can only say that |v | 2 is invariant. (The conservation of elec- 
tric energy — the invariance of s (TQ = (l/2) C | V | 2 , where then = 
K(z) — would have to be indicated Sy a z-dependent solution to 
Eq. (122).) On the other hand, an invariant form(l/2)Re VI* implies 
power conservation in any case, for the expression for the power 
does not depend on z-dependent parameters of the system. 

By introducing a still different set of K(z), we can parallel 
for nonuniform systems the technique discussecT in Section IV for 
determining the range of operators with specified appropriate 
metrics. But first let us point out that the results given there 
are not all restricted to uniform systems. If constant K.q and JK 2 
are specified, it is not in fact necessary to choose solutions .H 
to Eq. (Ill) to be constant. The arbitrary constants in the gen- 
eral Hermitian solution can as well be arbitrary functions of z. 
This leads by Eq. (110) to the complete range of systems, uniform 
and nonuniform, with operator I? having constant and K 2 as 
appropriate metrics. It is not necessary that there be any more 
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independent, constant metrics appropriate to an in this range, 
however, if R is not constant — unless there are more in the set 
K, J m- 1 implied by K, and K„ . 

More generally, if the conjugacy properties of JM(z,z ^) , the 
dimensionality of the eigensubspaces , and their organization into 
chains of generalized eigenvectors do not depend on z (except at 
discrete values) there will be p linearly independent solutions 
to the equation 

M + (z,z, )K(z)M(z,z 1 ) = K ( z ) (127) 

/W\ _L Aw Aw X ANA 

where p depends on the degeneracy in M(z,zi) exactly as was deter- 
mined for uniform systems. Note that for any lossless, nonuni- 
form system, with at least the one invariant s (JC) , K constant, re- 
presenting the power, there must by our earlier results be at 
least n independent, nonsingular solutions to Eq. (127) — for there 
exists an appropriate nonsingular K at every value of z. 


Equation (127) implies that 

x^ (z) K (z) x (z) = x^ (z. ) K (z) x (z, ) (128) 

/VVN Ana AV. AW X ✓*** tAW X 

For a given z, K(z) can be considered a constant metric. The s (JK) 
formed from solutions to Eq . (127) will then represent properties 

of the system that are invariant when transformed from zi to z, 
but the invariant properties will in general depend on z. But 
this suggests that we can define a nonuniform system in terms of 
its z-dependent invariant properties in a manner parallel to the 
way in which we defined a uniform system in terms of its z-indepen- 
dent invariant properties. 

Given two independent, nonsingular, Hermetian metrics Ki(z) 
and K 2 (z) , we can consider an M(z,z]_) for which these are appro- 
priate as represented by an equivalent uniform line from zi to z. 
The system operator for this uniform line can be denoted by Rq(z), 
with z as a parameter. Defining 

J (z) = K" 1 (z) K_ ( z ) (129) 

we obtain the Hermitian solutions to 

J + (z)H(z) = H (z) J (z) (130) 

MA AW AAa AW 
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Since H = & 1'&2 are nonsingular solutions to Eq. (130) , there must 
be at Teast n independent solutions. There may not be the same num- 
ber at every value of z, and the H.(z) may be otherwise not well be- 
haved, if the conjugacy properties of the eigenvalues or the de- 
generacy in J changes with z. But we can always choose at least 
n independent 1 solutions. We thus determine a range of Rq (z) by 

R n (z) = k: 1 (z) H ( z ) (131) 

amU /w\ 


Since Rq (z) is the operator for a uniform line, it is straight- 
forward to determine its matricant as the solution to the equa- 
tion 


dM„ (z' ,z, ) 

3 —as' '■■■ - 5o (z) “o (z ' ^i 1 * So 12 !' 2 ! 1 = l (132 > 


where Mq(z',zi) contains z as a parameter. Corresponding to the 
range "of JIq(z) is a range of Mq(z',zi) having J^i(z) andjj^^ 2 ) as 
appropriate metrics, and in particular this is true at z' = z. 
Thus Ki(z) and K 2 (z) are solutions to Eq. (127), with 


a <z - z i> -Mo lz '' z i > 


(133) 


and the M(z,z]_) are lumped representations for the range of sys- 
tems defined by these conservation laws. If the H(z) are suffi- 
ciently well behaved that the M(z,zi) are differentiable, then 
the latter are matricants from which a range of system operators 
can be determined by 

R(z) = M _1 (134) 
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This approach to system definition, in terms of metrics de- 
fined by Eq. (127), has the same limitations as the parallel 
approach for uniform systems. For this purpose, we shall find 
that the metrics defined by Eq. (122) are more significant. These 
are quite different definitions of K.(z) , for with JR. given by 
Eq. (134) , we find that solutions to 1 Eq. (127) must satisfy 


R + K - KR 

Ha* fw* 
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(135) 


and so do not in general satisfy Eq. (122). Incidentally, Eq. 

(135) indicates that we cannot obtain or characterize the solutions 
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to Eq . (127) by a parallel treatment of R, as we could for uni- 

form systems. The solution of Eq. (135) for requires the deter- 
mination of M, and Eq. (135) thus offers no real alternative to 
Eq. (127) in the study of nonuniform systems with this interpre- 
tation of the conservation laws . 


VI. USE OF DIRECT PRODUCT TO OBTAIN THE CONSERVATION LAWS 


The conservation laws governing a distributed, uniform or 
nonuniform system can be determined from another property of the 
Kronecker product and sum, if the matricant is known. Like the 
procedure described in Section III, and unlike that of Section II, 
this procedure does not require determining the eigenvalues, eigen- 
vectors, and generalized eigenvectors of M.(z,zq), or of R obtained 
by Eq. (46) . If only R is known, it is still sometimes easier to 
determine the matricant (e.g., by the method of projectors (ref. 3), 
for R constant) than the eigenvalues and modal matrix. If R is 
not constant, but satisfies the Lipschitz condition, the matri- 
cant can be determined at least approximately (refs. 3 and 11); and 
in either case this procedure offers a method of solution for 
Eq. (122), providing additional information beyond that given by the 
constant K. 

It is well known that the direct product of the matricants 
M A (z,z]_) and MB(z,zq) of two nxn operators A and B is the matri- 
cant of the Kronecker sum A (A,B) . If, that^is, *** 

am am am 


- 5 -“= — j AM , M ( Z-, , z, ) = I 
dz -~A 11 M 


dM 


df*= -3 Bm B ' J£ b (z 1 ,z 1 ) = £ 


(136) 


then, writing M^z^) x Mgtz^) as M A x M^z^), 
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— dz- = -Wt'2) Mb 1 ' Mb' 2 !' 2 !) = L (137) 
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It is easily verified that 


d (M A x M B ) dM A 
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Writing ^Jr=IMb in the first term and in the second term, 

and with the use of Eq. (82) and the definition, Eq. (84) , of 
A. (A, B) , Eq. (137) follows directly. Since the eigenvalues of 
Ma%B are products of the eigenvalues of Ma and the eigenvalues 
of Mb, and Ma and Mb are nonsingular, Ma x MeT is also nonsingular — 
as required 5? a matricant. The boundary condition follows from 
those in Eq. (136) and the definition of the direct product. 

We found in Section V that Eq. (122) can be put in the form 
of Eq. (125) , and that the same condition that ensures the exis- 
tence of a matricant for .R ensures the existence of n^ linearly 
independent solutions Jc" to Eq. (125). If any set of linearly 
independent solutions is arranged as the columns of an n2xn2 
matrix K" (z) , then K" (z) satisfies 


^ = -ii<s + '-s T ) r 


AM AV» 


(138) 


By Eq . (137), Eq . (138) is satisfied by 


K " ( z , z , ) = M , x m K" (z. , z. ) = I 

X T ^ T AM X X AM 

K “K 

AM AA 4 

,+ 


(139) 


where M , (z,z,) is the matricant of R and M _(z,z.) is the 

\ X flM A** X X 

I* T a* 

matricant of -R . To relate M , to the matricant M(z,z.) of R, 

aMA AVI X 

we take the transpose conjugate 4 of both sides of Eq. (39) and 

t -1 

postmultiply by (M ) , to obtain 


dM."^" ,„t.-i . u t D t (u tv ^ 

-g±- (M ) = jM R (M ) 

Cl Z A** /w\ a/v\ am 


(140) 


But 


§£ (M + ) _1 


-1 


= -M 
✓w» 


t d (M 1 ) 


dz 


(141) 


since (d/dz) [M^ (M^ ) = dl/dz = jO. Substituting Eq. 

Eq. (140) ancl premul tip lying by ^'(M+) - l, we obtain 


(141) in 


— = -jR + (M + ) 


-1 


(142) 
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is the matricant 


I. Thus (M^) ^ 

rtM HM 


(143) 


Similarly, we find that 

M = (M T ) _1 (144) 

fW* 

Thus 

K" = (M + ) _1 x (M T ) -1 (145) 

rtrtA /Wt A* 

A nonsingular solution to Eq. (138) is obtained from the 
matricant of .R by Eq. (145) . The columns of K" thus form a set 
of n^ linearly independent solutions kV to EqT (125) . 

Any linear combination ^afk^ of these solutions, where the 
ai are independent of z, is also a solution. An elementary trans- 
formation, as defined in Section III under "Testing the Kronecker 
Sum," on the c olumns of K" , will result in a new matrix with 
columns that are linear ’combinations of the columns of K" . Further- 
more, it is a property of elementary transformations that they do 
not change the rank of a matrix. The columns of the new matrix, 
therefore, will also form a set of linearly independent solutions 
to Eq . (125) . 

All constant solutions to Eq. (125) will be linear combina- 
tions^of the columns of K " , obtained by taking the direct product 
of (M t ) - 1 and (M T )~1. The constant solutions can be obtained by 
successive elementary transformations on the columns of K" , with 
the assurance that after each step we will be left with a set of 
linearly independent solutions . A scalar multiplier employed in 
an elementary transformation must be chosen independent of z, and 
the functional forms of the elements of Kj' will serve as a guide 
in making elementary transformations that will display the constant 
.k" . If the required transformations are not evident, we can sys- 
tematically obtain as many constant elements as possible in the 
first row. The constant k" must then be linear combinations of 
the columns headed by the constant elements, and by elementary 
transformations on these we can obtain as many constant elements 
as possible in the second row. Continuing in this way, we will 
obtain the largest set of linearly independent k" that can be 
chosen constant. 
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The constant metrics K appropriate to R are obtained by re- 
arranging the elements of the constant k." as n*n matrices in 
accordance with the convention we have adopted. As the constant 
solutions to Eq. (122) , they will be the constant solutions to 
Eq. (8). The remaining independent solutions J^(z) to Eq . (122) 

are obtained similarly from the remaining k" (z) , and for nonuni- 
form systems these may imply additional conservation laws. Note 
that if K is a solution to Eq. (122) , K+ is also; and again the 
solutions can be chosen Hermetian. 

Example 1 ; 


The operator R for a uniform transmission line, given by 
Eq. (3) , has the matricant 


M (z , 0) 

M A 




By Eq. (145) 


we obtain 


(146) 


K" 



Bz cos 
B z cos 


Bz 


Bz 



j =i- sin Bz cos Bz 

z o 

2 o 

cos B z 
2 

sin B 2 

-jZq sin Bz cos Bz 


-j sin Bz cos Bz 

Z 0 

■ 2 u 
s l n B z 

2 

COS I* z 

j Z sin | z cos B z 


1 . 2 . 

— 5- sin Bz 

Zo 

-j ^ sin Bz cos 

Z 0 

j sin Bz cos 

z o 

2 

cos B z 



Writing the i ' th column 



(147) 


(148) 
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we see that (l/M^k^+k^ ) gives the elements of K 1 , and (l/4Z 0 )(k^+ZQ^^ ) 
the elements of K 2 , in Eq. (60). If k.£ and Jk" are replaced res- 
pectively by (l/4)(k 2 +k 2 ) and (l/4Z Q )(k^+z2k^ ) , the result of elemen- 
tary transformations, it is clear that we cannot then obtain 
another column with constant elements by replacing or k.^ by a 

linear combination a-.klj+a.k" , a 0 and a. constant. 

3»*3 4-^4 3 4 

Example 2 : 

The operator R for an exponentially tapered transmission line 
is of the form ^ 



(149) 


where L=L 0 e az is the inductance per unit length, C=Coe is the 
capacitance per unit length, Zq = ^ Lq/Cq , and 3 = co ^ LC = w V L 0 C 0 • 


-az . 


The matricant is 


/(l/2>; 


M ( z , 0 ) = 
<w> 


( 


1 a 


cos yz- j - sin yz j -j 


)- 


, 3Z 0 (l/2>z 


Y 


sin yz 


V 


-j 775 — S 1/2 ) JZ Sin yz 




yZ, 


(■ 


cos yz+ — sin y: 
z y 

(150) 


2 2 12 
where y =$ - j a . 

vectors ; 


By Eq . (145) we obtain K" , with the column 

I'M 
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Clearly (l/4) (k^J+k^ ) gives the elements of 


K n 




( 152 ) 


and so the power is conserved (see Eq. (6)). Successive elementary 
transformations, using the systematic procedure we have described, 
reveal that there are no constant metrics independent of . But 
if we substitute e az = (l/3Zo)r]_2 and e~° z = (Zo/3)r2i where these 
appear explicitly, and in the same manner seek linear combinations 
of the that are not explicit functions of z, we find that 

(B/Zq)^ + jcrk^ + $ z 0%;'£ gives the elements of ^ 


K 

n/V> 




(153) 


from which we can 


form the Hermetian solution to Eq. 


( 122 ) , 


£ 2 


_1 

4oo 


(K+K + ) 




(154) 


Equation (154) implies the invariance of 


s (K 0 ) = x K 0 x = 

AV'Z /W\ -V** 


C | V | 2 + | L | I | 2 


+ j j £ (V*I-VI*) (155) 


Thus the total stored energy plus a/to times the imaginary part of 
the complex power is conserved. 

If k-7 and k '2 are replaced respectively by (l/4 ) (k^+k^ ) and 
(3/Z^)kj + jok^ + 3 Zgk^ , we are left with k^ and k£ . Two additional 
independent conservation laws can be derived from these on sub- 
stituting z = (l/a) In (ri2/8Zo) or z = (l/a) In (B/Z q^i) everywhere ; 
but they do not have an obvious physical interpretation. 
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VII. GROUP REPRESENTATION WITH THE SYSTEM INVARIANTS 


We observed in Section IV that the constant, nonsingular me- 
trics appropriate to a system operator define this operator only 
within a range. It is also known that the matricant has group 
properties , and the group representation it provides is not unique 
(refs. 3 and 7). Rather than considering, in the usual manner, 
group properties of the matricant as defined by appropriate, non- 
singular, constant K, we shall find that the constant K, singular 
and nonsingular, are 1 embodied explicitly in another group represen- 
tation. That they are neither necessary nor sufficient for such 
a representation, however, will indicate the extent of their con- 
tribution in defining a system. We shall see what is needed to 
complete the representation, and shall find that it has both phy- 
sical interest and some further group properties. 

A nonempty set S is said to be a ring under a Lie product 
if it is closed under defined binary operations of addition (+) 
and multiplication (•)/ and provided, for arbitrary elements si, 
s 2 ' s 3 £ S : 

1. ( s i +s 2^ + s^ = + (S2+S3) . 

2. s ± + s 2 = s 2 + Sl . 

3. There exists an additive zero element sqsS such that, for 
each s^eS, s^ + Sq = s^. 

4. For each s.eS there exists -s . eS such that s. + (-s.) = 

1 1 11 

s 0 * 

5. s^ • (s 2 *S 2 ) + s 2 * (s^'s^) + * ( s ^ • s 2 ) = Sq (Jacobi 

identity) . 

<5. S;L • (s 2 +s 3 ) = S;L * s 2 + s 1 • s 3 . 

7. s 1 * s 2 = - (s 2 * s 1 ) . 

A multiplicative operation having properties 5, 6, and 7 defines 
a Lie product; and an example is given by the commutator of ma- 
trices A and B, defined by 

A • B = [A , B] = AB - BA (156) 

iW\ rw> M* jV* AW 

where AB^ is interpreted as ordinary matrix multiplication. 
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Let S be a subspace of a linear vector space over a field F 
such that S is closed under addition over V, where V is a sub- 
field of F. (If SjL»S 2 eS and vi,V 2 £V, then viS]_+v 2 S 2 eS . ) Then if 
S is a ring under a Lie product, S is said to be a Lie algebra 
over V. 

Let K be a given constant metric, and let be an operator 
defined over the field F of complex numbers. Let (+) denote 
ordinary matrix addition, and the commutator be taken as the Lie 
product. Then it is easily verified that S = -jR, where R is the 
range of operators for which mi 


R! K = KR . , 


R. eR 

VWil ^ 


( 8 ) 


is a Lie algebra over the field V of real numbers. The metric 
thus defines a Lie algebra for S. Furthermore, the specification 
of additional independent, constant metrics will define a Lie sub- 
algebra . 

A nonempty set G, together with a binary operation («), is 
called a group provided: 

1. If g 1 ,g 2 eG, then q^oq^eG (the group property) . 

2. ( gi *g 2 ) o g 3 = g 1 o (g 2 »g 3 ) , for all g 1 ,g 2 ,g 3 eG 
(associativity) . 

3. There exists an identity element leG such that, for all 
g^c, 

g i _« l = i • g ± = g ± 

4. For every g^eG, there exists an inverse g^^sG, such that 

-1 -1 , 
g . o g . = g . « g . = 1 


It can be shown that the identity and the inverse must be unique. 
If G has only the group property, G is called a groupoid. 


Again let K be a given constant metric. Then it is easily 
verified that M(z,z.), where M is the range of matricants for 

. _L avv 


which 


mTkm. 


K , M . eM 

AM /ml AW 


(36) 
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is a group under ordinary matrix multiplication. It is a contin- 
uous, finite group, in that it can be parameterized in terms of 
continuous values of a finite number of parameters. We have seen 
that if Mi is a matricant of , then Eqs . (8) and (36) have the 

same constant solutions K. The metric K thus defines a group M, 
and an associated Lie algebra for _S. The group M is known as a 
Lie group. The specification of additional independent, constant 
metrics will define a Lie subgroup and an associated Lie sub- 
algebra. (A rigorous treatment of what is meant by continuous, 
finite, and Lie groups is beyond the scope of this discussion. 

There are many excellent references; e.g., refs. 12 and 13.) 

The group M can also be thought of as a matrix representation 
of an abstract group G of transformations characterizing a range 
of systems. A homomorphism of a group G into a group or groupoid 
H is defined as a mapping g -*■ ga of each geG to h = ga, heH, satis- 
fying the condition 


(g 1 °g 2 ) a = (g^) ° (g 2 °0 * g 1 »g 2 eG 


(157) 


The binary operations for G and H need not be the same. We de- 
note by Ga the set of ga for all geG. Then Ga is a subset of H, 
and the mapping is onto Ga (i.e., for each heGa, there exists a 
geG such that ga = h) . It can be shown (ref. 14) that if a is a 
homomorphism of a group G into a group or groupoid H, then Ga is 
a group. The group Ga (or the mapping) is called a representation 
of the group G. If the mapping is one-to-one (i.e., gqa = g£cx 
implies gi = g2), the mapping is called an isomorphism, G and Ga 
are said to be isomorphic, and as the groups are essentially the 
same, we write G - Ga. The representation is then said to be 
faithful. Thus G = M can be interpreted as a statement that M is 
a faithful matrix representation of an abstract group G. ^ 


Let M,' be the set of all nonsingular, differentiable, nxn 
matrices Mj_(z,z ]_) , where Mj_ (zi , Z]_) = jl. Again it is easily veri- 
fied that these form a group under matrix multiplication. Then 
G' = M' is the abstract group encompassing all systems having n 
degre’es of freedom that can be described by a matricant, M' = G'a 
is a faithful representation of this group, and the set wTth ele- 
ments — j Rj_ = (dM-j_/dz)Mjl , MieM', is the Lie algebra associated with 
M ' . ‘ 


Assume P* = G’6 and Q' 
sentations of G ' , where for 


= G'y are two faithful matrix repre- 
all g^eG ' , 


g. 8 = P. , P . eP ' 
g, Y = Q. , Qj eQ ' 

1 M 1 I** 


( 158 ) 
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and the P^-eP' and Q- eQ' 1 are of order n*n. Let N' = G'n be defined 
by the mapping, for all g^eG' , 


9i n = ( gi 6) x (g^) = Mi * ii = Si' 


(159) 


where £i x Qi is the matrix direct product defined by Eq . (77) . By 

Eq. (82) ^and since 8 and y are homomorphisms , we can write 

( g i n ) • <* 2 n> = ( Hl x Sl ) ( ^2 X 22 ) = <SlS2> x ( &l£2> 

= [ (g-j_3 ) ° (g 2 3) ] * [ (g 1 y) « (g 2 Y) ] =[(g 1 °g 2 )3] 


[(g 1 °g 2 )Yl = (ii'-g^n 


1 y 2 ' 


for arbitrary gp,g 2 eG'. Thus N ' is a groupoid under matrix multi- 
plication, and n is a homomorphism of the group G 1 onto N' . It 
follows that N' is a group and a representation of G' — but not 
that the representation is faithful. 

Let G = M be the group defined by some set of constant me- 
trics K. G is a subgroup of G' , and M = Ga is a subgroup of M' . 
Then P*"= G8 and Q = Gy , defined by Eq . (158) for gieG, are sub- 

groups of P/ and**Q ' respectively, and faithful representations of 
G. This follows ‘since 8 and y are homomorphisms of the group G 
into the groups P' and Q', and since the representations P' and 
Q' are faithful."** Similarly, n is a homomorphism of G into N'; 
and N = Gr) , defined by Eq. (159) for gj_eG, is a subgroup of"*N' 
and a representation (not necessarily faithful) of G. 

Now we note that taking the inverse of the transpose conju- 
gate of all members of the group M.' , or the inverse of the trans- 
pose, constitutes a one-to-one mapping of M' onto itself. Further- 
more, each of these mappings is a homomorphism, since for arbitr- 
ar Y M-l, M 2 eM', 

GW 1 = MT 1 

(el )' 1 ^)" 1 - [MV 

If, that is, we define P' = M'p and Q' 

/W\ aao /\aa 


= M't, where for all M.eM', 

ftAA AM 
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( 160 ) 


/ 

+ \-i 


M . p = ( IN 

( 

k) 

t \ -1 

= P, , P. eP' 

/Wsl AM 

M. T = I IS 

X \a 

£ 

= Q-: , Q.eQ' 

AmJ- Mi AAA 


then p and t are isomorphisms. Accordingly, M' =" P 1 and M' = Q'. 
We can also write *** *** 


P' = G'ap = G'P 

/MA 

Q 1 = G ' ax = G ' y 


(161) 


and since a is an isomorphism as well, 8 = ap and y =ax are iso- 
morphisms. Thus G" *= j?" and G' = Q', and P/ and <^' are faithful 
representations of G 1 . Also, thenT J? = Mp and Q = fix, defined by 
Eq. (160) for M-j_eM, are faithful representations' of G. It follows 
that N' = G ' ri ^where for all g . eG ' 

r i 

g.ri = P. x Q. = (m?\ x = N. , N.eN' (162) 

y x ' «*x U-iy V* 1 / rt * 

is a representation of G ' , and the subgroup N = Gri is a represen- 
tation of G. ^ 

Thus the K" of Eq. (145), deriving from all M^eMj , are seen 

to form a group N' , and the .K" deriving from M-leM a subgroup 14, 

that are moreover 1 representations respectively of G' and G. But 
are these representations faithful? Inversion of Eq. (138), 



(K") 


AAA 


-1 


t T 
A (R,-R A ) 


(163) 


t T 

shows that K" determines the Kronecker sum AJR^ , - R, ) . Inspection 
of Eq. (97 ^reveals further that if A (RT , -R+) A(R~, -rT) , then 

r 2 = j^l + a i,f a real. Thus, while strictly the representations 
are not faithful, K" determines R within the identity. 

O 

The columns k" of a given K"eN' or K"eN are generators of n 
independent invariants s (K) that define the system operator with- 
in the identity. Not only do a set of constant K define a group, 
and a range M of that constitutes a matrix representation 

of the group7 but** written as Jc" they are explicitly displayed as 
linear combinations of the columns in another matrix representation 
of the group. The system is not defined by the constant K alone. 
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and we see that the IC ( z ) are essential to complete the group- 
theoretical representation. Moreover, this extension of the set 
of constant is not just formalistic; for we have seen that, for 
nonuniform systems in general, invariants obtained from the K(z) 
reasonably admit of interpretation as conservation laws. From 
this point of view, the class of uniform systems can be considered 
as an atypical case, represented by a "singular" surface of N' — 
or if in the range defined by the constant K > , by a singular 'sur- 
face of N. This singularity represents an inability to interpret 
the invariants as conservation laws, but not an inability to 
obtain pre-images from the representation. 

If K"eN (or K"eN') is the matricant of A(R^,-R T ), then -jA(R+, 
-R ) is a member of the Lie algebra associated with N (or N 1 ) . 

This can be verified directly. To verify closure under commuta- 
tion, let -jR 3 be the commutator of two operators -jRg and -jR 2 
in a Lie algebra defined by a set of constant K, 


= [ "=>Si'-jS 2 ] 


(164) 


Then -jR 2 is also in the algebra. It is easily checked that 


' j A (R* 



-jA (R^,-R^) ] 

AA A*v Z. mZ 


- jA (R^t 



(165) 


4* rp 

The columns of the matricant of the operator A(Ro,-Ro), therefore, 

+ __ ANA 

T 4* m 

witn -]ji (Ki , -R 3 ) obtained by commutation of - j A (R-i , -R n ) and 

, XVAA. /W.-L- «**-L 

*f* T 

“ j 4 ( R ? / -R 2 ) ' will also yield in linear combination the k" formed 
from the specified K. The other conditions for the K to define a 
Lie algebra for -jA(R^,-R T ), obtained by Eq. (163) for K"eN, can 
likewise be verified? ^ ** 

Kronecker sums in this particular form, then, might be des- 
cribed as providing a representation of conservation laws onto 
their algebra. 

If the constant K define a group, what can be said of the 
range of systems having an invariant s (JK) , where .K = K(z)? The 
range of matricants will be the solutions to Eq. Tl24?, 

M + K(z)M = K ( Z-. ) (124) 

Am /Wi aw v) X 


with the boundary condition M(zg,zg) = I. If we define a binary 
operation for arbitrary in this "range by 
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(166) 


M 0 = M, 

/WO ^M_L 


M„ 


M„M, 1 M„ 
a«2 


it is readily verified that M 3 also satisfies Eq. (124). Thus 
the matricants in this range form a groupoid under the defined 
binary operation. The binary operation is not associative, how- 
ever, and the groupoid does not contain an identity element. Thus 
it is not a group. The mapping -*■ (Mj ) - 1 x (mT )~1 is to a 
groupoid under this same binary operation. Again it is a homomor- 
phism, and the form a representation of the groupoid defined 
by the invariant s (K) . 


By way of example, consider the exponential line of Example 2 
in Section VI. The metrics K 3 and ^(z) °f Ec 3 s • (152) and (154) 

are solutions respectively to Eqs . (36) and (124), for the matri- 

cant of Eq. (150). Let us denote this matricant by M]_. But K 3 
and £2 are also solutions for ^ 



Thus and are in the group defined by JK 3 , and in the group- 
oid defined by K 2 . We can obtain another member of the group. 


az 


( c 


1 a . 


e (cos yz - ^ — sm yz] -j — sin yz 


■) 


ez 


0 , 


M-j = M,M 0 = 


YZ, 


sm yz 


-az 




,1a . 

e \cos yz + 2 — sm yz 


z 

/ 


(168) 

t 

which satisfies M 3 K 3 M 3 = K 3 , but which is not a member of the 
groupoid (M^K.2 ( z );^3 ¥ K2(0))* We can obtain another member of 
the groupoid. 
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/, /0 n / . 1 a . \ . BZ 0 (l/2)o z . 

e (l/2)ozlcos yz + j y sin Y z ) 3 "tjp - e v ^ si 


m yz 


M /t =M 0 M 1 1 M 0 = 

AM ft **X W 


. 3 -(1/2)0 z . 

jdr- e ' ' ^ sin yz 


j"0-/?^ z /^. os y Z - ik £. sin yz 

V 2 Y / 


(169) 

which satisfies M 4 & 2 (z)M 4 = K 2 ( 0 ), and it is seen that M 4 is also 
a member of the group. ~ 

Members of a groupoid, defined by some £ 2 (z), which are gen- 
erated by other members Mj and M 2 that are also members of a group 
defined by constant K lf will always be members of the group. For 
if MjKiMi = Ki and mJk;|M 2 = Ki , then 


(M, °M 0 ) + K, (M, °M 0 ) = M** w 

(Wl Av\Z /+Al_ A*X yw>Z 


= M 


ifi 1 ) 

(“I 1 ) T £12/; 


K-| M„M-, 1 M 0 
**X **Z »*X a*v- 


M 0 = MiK,M 0 = K, 

•maZ aw* z *aX 


(170) 


The role of the conservation laws in defining a linear sys- 
tem is thus clear. A system having n degrees of freedom is de- 
fined by its n2 independent invariants s ( K ) . The constant K, if 
any, including singular obtained from linear combinations of 
the columns of the nonsingular matrix K/ , define a group — the 
intersection of the groups defined by the constant separately. 
The £(z) , which we have seen in general lead to valid conserva- 
tion^laws, define a groupoid — the intersection of the groupoids 
defined by the I^(z) separately. By Eq. (170), intersections of 
groups and groupoids define subgroupoids . The intersection of 
all groups and groupoids is a subgroupoid that defines the sys- 
tem operator within the identity. 

The n conservation laws cannot, however, be specified arbit- 
rarily. A specified K" must be chosen in the group N' , which is 
a subgroup of the group formed by taking all possible matrix 
direct products Pi x Qj of the J?i e L' and Qi £ Q' * ( The latter is 
a representation“of a^group- theoretical direct product, a slightly 
different concept that we have not discussed.) This in turn is a 
subgroup of the group of n2 x n2 matricants. That is to say, the 
groups and groupoids defined by some sets of are disjoint, and 
others are included in one another. 
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The M-^eM, the group defined by the specified constant K, 
evaluated at zj_ = 0 and any z = L, is also a group. The mapping 
is a homomorphism of the group with elements Mi(z,0), M^eM, onto 
the group with elements Mi(L,0), for all M.j_(z,0)eM. It is not 
an isomorphism, since for any Mi(L,0), M-^ (z , 0 ) eM^there will exist 
Mj(z,0)eM, Mj ( z , 0 ) ^ Mi (z ,0) , 'such that^M j (L , 0 )^*= Mj_ (L , 0 ) - the 
equivalent distributed system for a lumped system is not unique. 
Also the mapping depends on the choice of z = L. But for each 
z = L ^ 0, the Mi(L,0), Mi(z,0)eM, will include all nonsingular 
matrices K-Unitary with 'respect to the specified K, and will pro- 
vide an unfaithful representation of G . It is thus not necessary 
to determine Mi(z,0)^M that map to the group of ^nonsingular 
Unitary matrices — such pre-images as satisfy Mj_ (z , 0 ) K (z)M^ (z , 0 ) = 

K (0) , with K (L) = K ( 0 ) , and where the specified K are^among the 
K ( 0 ) . FurtKermoreT for each Mj_(L,0), (z ,0) eM, "''among the Mj(z,0) 
eM that map to Mi(L,0) there will exist matricants of constant 
operators R. It is thus not necessary either to determine all 
MieM — which requires obtaining the matricants of all R, including 
R(z) , with -jR in the associated Lie algebra — in order to obtain 
the group of lumped-system operators Mj_(L,0), M-j_ ( z , 0 ) eM , defined 
by the specified constant K. ^ 


VIII. CONCLUSIONS 

A distributed system described by a matricant M(z,z]_) that 
is a member of various groups and groupoids , with an operator -jR 
in corresponding algebras, can be described also by a matricant 
K"(z,zi) with the same group properties, and an operator -jA(R+, 
~RT) in corresponding algebras . The latter representation ex- 
hibits explicitly the metrics that define the conservation laws 
governing the system. It establishes and elucidates the signifi- 
cance of all the system invariants in characterizing a system. 
Their interpretation as conservation laws can be supported also 
by physical arguments. 

This representation can be viewed as the end result of a 
theory developed along parallel lines to existing theory for lin- 
ear, uniform, homogeneous systems. It provides new insights and 
techniques for the analysis and synthesis of uniform systems in 
terms of their conservation laws, and extends to a larger class 
of systems . 

The representation provides a model that can be studied for 
synthesis from any number of specified metrics. Synthesis may in 
principle be carried out directly by filling in the columns of 
Kj'(z,zi), since this is an equivalent matricant for the system. 

The problem is that this must be done under the constraint JK^'eN' . 
But an intriguing feature of this approach is that one is specify- 
ing the system and its conservation laws simultaneously and 
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explicitly in forming the matricant K"(z,z^). The constraint 
under which this must be done is in fact a fundamental property 
of the group that is revealed by this representation. There are 
inherent restrictions on the conservation laws that can be obeyed 
together by a system. 

Mode coupling in systems is a particularly appropriate area 
of application for this theoretical model. Considerable physical 
insight can be obtained by focusing attention on conservative ex- 
changes other than the power exchange between modes. As should 
be expected, the physical information to be specified or obtained 
is increased if the coupling or the coupled systems are nonuniform, 
although the number of linearly independent invariants from con- 
stant metrics will in general be fewer. The representation pro- 
vides a basis for characterizing, classifying, and synthesizing 
such systems. 
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